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Volume IX (1958) de ZAMP paraitra en deux parties, 
VOL. IXa 


comme suit: fasc. 7 le 25 janvier, fasc. 2 le 25 juillet, fasc. 3 le 25 septembre et 


fasc. 4 le 25 novembre 1958, tandis que le 
VOL. IXb 
FESTSCHRIFT JAKOB ACKERET 


sera un fascicule double No. 5/6, de 777 pages et paraitra au mois de juin 1958. 
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Kennzeichen atmosphiarischer Eispartikeln 
2, LEW 


Hagelkorner 


Von Roranp List, Weissfluhjoch-Davos?) 


1. Einleitung 


In die Augen fallend waren seit jeher die vielfaltigen, zum Teil auch phan- 

tastischen Formen und Grdéssen der Hagelkérner. W. TRABERT [1]?), der 1899 
eine vorziigliche Zusammenstellung der damaligen Kenntnisse tiber die Bildung 
_ des Hagels verfasste, erwahnt zum Beispiel Korner mit aufgebauten, bis 14 mm 
langen Kristalliten; C. KASSNER [2] beschaftigt sich mit der Frage des gréssten 
Hagelkornes der Welt, und S. Ers [3] beschreibt eine Schlosse mit 810 cm* 
_ Verdrangung. 
: Die Erkenntnis, dass ein schalenartiger Aufbau von den in das Eis einge- 
bauten Luftblasen herriihrt, ist bis auf WALLER und HaArtIne [4] zuriickzu- 
verfolgen. Diese Autoren konnten auch feststellen, dass die Schalen der Hagel- 
k6rner «aus einer Anzahl kleinerer, aneinander gereihter Theilchen, kleine Eis- 
_ pallchen oder Eiszellen, die man deutlich unter dem Mikroskop beim Schmelzen 
wahrnimmt», bestehen. Diese Beobachtung wurde hundert Jahre spater durch 
R. List und M. pe QuERvatn [5] bestatigt, da sich bei Betrachtung von Diinn- 
schnitten in polarisiertem Licht diese «Eisballchen oder Eiszellen» mit Eisein- 
kristallen identifizieren liessen. Die damals erstmals auf Hagelkérner ange- 
wendete Diinnschnittechnik [6] wirkte befruchtend, indem die damit erreich- 
baren Resultate verschiedene Vermutungen bestatigen, in anderen Fallen aber 
neue Ergebnisse zeitigen. 

Die folgenden Ausfithrungen basieren auf Sammelproben von iiber 40 Hagel- 
fallen der Jahre 1953 bis 1957. , 


1) Eidgendssisches Institut fiir Schnee- und Lawinenforschung. 
2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis Seite 234. 


218 Ro.anp List ZAMP 


2. Charakterisierung 


Aussehen. Hagelkorner sind teilweise oder ganz opake Eispartikeln mit einem 
Durchmesser von rund 5 bis 150 mm und mehr. Die auftretenden Formen sind 
dusserst vielfaltig; der Aufbau ist 6fters schalenartig mit einem Wechsel klarer 
und opaker, weisser Schichten. 

Entstehung. Hagelkérner entstehen aus Frostgraupeln oder Klareispartikeln 
durch Anlagerung und Gefrieren von angeschwebten Wassertropfen [7] (vgl. 
auch mit MAson [8}). ; 

(Die in der Meteorologie vielfach empfundene Unsicherheit in der Unter- 
scheidung von Frostgraupeln und Hagelkérnern einzig nach einem Grenzdurch- 
messer von 5 mm [9] kann durch Beachtung einer prinzipiellen Differenz be- 
hoben werden: Frostgraupeln sind Partikeln mit einer zonenweisen Verglasung 
und stellen Verdichtungsformen verschiedener Stadien dar [7]. Sobald der 
erste auftretende Verdichtungsprozess, ausgehend von der Reifgraupel, be- 
endet ist, das heisst, sobald die angestrémten Wassertropfen ausschliesslich zur 
Volumenvergrésserung dienen und nicht mehr zur Auffiillung des «Reif»- 
geriistes, so ist bereits von Hagelkérnern zu sprechen.) 


3. Ausserer Aufbau von Hagelkérnern 
3.1 Ursprung 


In Bestatigung des im Artikel Graupeln als Wachstumszentren von Hagel- 
kornern |7| erwahnten Schemas iiber die Entstehung und den Ursprung von 
Hagelkérnern, sei auf die Figur 1 verwiesen, die im Innern eines rund 5 cm 
grossen, ellipsoidischen, in Richtung der kleinsten Achse beidseitig eingedriick- 
ten Hagelkornes (Nr. 57.7) deutlich eine urspriingliche, kugelsektorartige Reif- 
graupel erkennen lasst. Eine Vergrésserung der Graupel zeigt Figur 2. 


Figur 1 Figur 2 
Hagelkorn 57.7 mit eingebauter konischer Kegeliges Zentrum des Hagelkornes 57.7. 
Reifgraupel. 
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3.2 Formen und Gréssen von Hagelkirnern 


Die ausserordentliche Vielfalt der beobachteten Formen lasst folgende 
_ Unterteilungen zu: 


Hauptform: 

a) Kugel; 

b) Ellipsoid: 
Rotationsellipsoid ; 
Ellipsoid, dreiachsig ; 
Ellipsoid, einseitig in der kleinsten Achse eingedriickt ; 
Ellipsoid, beidseitig in der kleinsten Achse eingedriickt ; 
runde Platte; 
runde Platte mit Wiilsten; 
ovale Platte; 
ovale Platte mit Wiilsten; 

c) Kegel, Kegel auf Kegel; 

d) Schnitz; 

e) unregelmassige Form; 

 f) + teilweise oder ganz ebenflachige K6rper; 

g) zusammengeballte Partikeln ; 

h) Splitterformen. 


Oberflache: 


aa) gleichmassig gekriimmt ; 

bb) unregelmassig ; 

cc) protuberanzenartig ; 

dd) stachelig; 

ee) mit aufgebauten Kristalliten. 


Zu den zusammengeballt aussehenden Kornern ist zu bemerken, dass der 
Diinnschnitt unter polarisiertem Licht das einzige sichere Kriterium ergibt, um 
zu unterscheiden, ob es sich um eine reelle Zusammenpappung verschiedener 
Eispartikeln handelt oder nicht. Alle scheinbar durch Koagulation entstandenen 
Korner, die auf Weissfluhjoch untersucht wurden, entpuppten sich als ganz 
gewohnliche Korner mit esnem Kern! 

Die in einem einzigen Hagelschlag auftretenden Variationen der Korn- 
formen kénnen dusserst vielfaltig sein, wie die Untersuchung der Sammelprobe 
vom 24. Juni 1953 zeigt [5], oder auch mehr oder weniger homogen gemass 
Figuren 3 und 4, die in erster Linie ellipsoidische Formen darstellen. 

Gréssenverteilungen von Hagelkérnern konnten nur zwei, die den gleichen 
Hagelschlag vom 5. August 1955 in Balsthal betreffen, gemessen werden; sie 
sind in Figur 5 aufgezeichnet. 
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Sammelprobe des Hagelschlages vom 21. Juni 1 
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3 4 5 6 i 8 9 10 "1 12mm 
Durchmesser 


Figur 5 


Korngréssenverteilung des Hagelschlages vom 5. August 1955 in Balsthal. Totalzahl 990 bzw. 
437 Partikeln. 


3.3 Dichten von Hagelkornern 


Samtliche iiber Auftrieb in Petrol oder Tetralin und Wagung bestimmten 
Dichten liegen zwischen 0,8 und 0,9 g/cm’. Die untere Grenze diirfte sich im 
allgemeinen noch nach kleineren Werten hin verschieben, da die gesammelten 
Schlossen vielfach vor der Einbringung in Kiihlschranke oder -réume ober- 
flachlich mehr oder weniger auftauten. 


4, Anordnung, Form und Grdésse eingebauter Luftblasen 
4.1 Schalenanordnung 


Schalenanordnungen von Luftblasen in Hagelkérnern sind Kennzeichen 
fiir eine Variation der Wachstumsbedingungen. Die Form der Schalenflache 
kann kugeligen (Figur 6), ellipsoidischen (Figur 7), kegeligen (Figur 8) oder un- 
regelmassigen Charakter (Figur 9) aufweisen. Die Grossen der schalenbildenden 
Blasen liegen im Mittel zwischen 2 und 100 p, ihre Form ist kugelig oder un- 
regelmassig, sofern sie teilweise vom urspriinglichen Korngeriist begrenzt 


werden. 
4.2 Rethenanordnung 
Vielfach sind in Hagelkérnern neben Schalenanordnungen oder auch allein 


Reihenanordnungen von Luftblasen festzustellen (Figur 10), deren Auftreten 
jedoch nicht abhangig ist von einer Anderung der Wachstumsbedingungen. 
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Figur 6 Figur 7 


Figur 8 Figur 9 


Blasenschale unregelmassiger Form, 
Kor 57°12 .Da- 


Figur 10 


Figur 11 


Reihenanordnung von Luftblasen, Luftkanale, Korn 57.16, D1. 
Korn 57.10, D 1. 
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_ Thre Richtung ergibt direkt diejenige des Wachstums. Die Durchmesser der 

_ reihenbildenden Blasen liegen meist innerhalb 2 und 200 u. Diese Gréssen wer- 
den aber iiberschritten, wenn es sich um radiale Kanidle handelt, die nach- 
_ traglich am Ende zugefroren sind (als Beispiel siche Figur 11). 


5. Kristallographischer Aufbau 
5.1 Bedeutung der Diinnschnittlage im Korn 


Bei der Auswertung von Einkristallanordnungen ist es von Bedeutung, aus 
welcher Lage der Diinnschnitt aus dem Hagelkorn herausgeségt wurde. Ver- 
gleichbare Verhaltnisse sind nur erhaltlich, wenn auf Grund der meist zuver- 
lassigen dusseren Symmetrie Hauptschnitte, das heisst Schnitte durch das 
Zentrum, gelegt werden. Kontrollschnitte helfen dabei diese Hauptschnitte zu 
bestatigen. Zur Demonstration der dabei auftretenden Eigenheiten seien ver- 
schiedene Diinnschnitte durch zwei Korner gleichen Typs dargestellt. Ihre 
Lagen innerhalb der beidseitig eingedriickten Ellipsoide seien anhand von 
Figur 12 aufgezeigt. 
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Figur 12 
Darstellungen der Schnittlagen durch die Korner 57.13 bzw. 57.10. 


Diinnschnitte von Korn 57.13 (Figuren 13 bis 18), die parallel zur kleinsten 
und gréssten Hauptachse ausgefiihrt wurden, geben indirekt fiir den schaligen 
Kornaufbau Hinweise; wahrenddem Schnitte durch Korn 57.10 parallel der 
mittleren und gréssten Hauptachse (Figuren 19-24) direkt den gleichen Schluss 
erlauben. Bei Korn 57.13 ist das Zentrum sicher kegelig, wahrenddem bei der 
andern Partikel eine kugelige Graupel als Ursprung angenommen werden kann. 
Zu Figuren 13 und 14 muss noch die Bemerkung gemacht werden, dass aus 
technischen Griinden der Diinnschnitt nur durch einen Teil des Kornes gemacht 
werden konnte; zur besseren Orientierung ist jedoch die Symmetrieachse ein- 
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Figur 13 Figur 14 
Korn 57.18, D 1. Korn 57.18, D 1, mit polarisiertem Licht. 


Figur 15 Figur 16 
Korn 57.13, D 2. Korn 57.13, D 2, mit polarisiertem Licht. 


Figur 17 Figur 18 


Korn 57.18, D 4. Korn 57.13, D 4, mit polarisiertem Licht. 
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Dinnschnitte von Hagelkérnern in polarisiertem Licht 


Dia Nr. 1, Korn 57.A. D1 (1:1,15). Dia Nr. 2, Korn 57.B, D 1 (1:1,15). 


Dia Nr. 6, Zentrum D 1, Korn 57.7 (2,1:1). 


Dia Nr. 5, Zentrum D 2, Korn 57.7 (2,1: il) 
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Diinnschnitte von Hagelk6érnern in polarisiertem Licht 


Dia Nr. 11, Zentr corn 57.11 (2,1: ia N 
ia Nr. 11, Zentrum D 1, Korn 57.11 (250 aL). Dia Nr. 12, Zentrum D 2, Korn 57.14 (2,1: i) 
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Figur 19 Figur 20 
Korn 57.10, D1. Korn 57.10, D 1, mit polarisiertem Licht. 


Figur 21 Figur 22 
Korn 57.10, D 2. Korn 57.10, D 2, mit polarisiertem Licht. 


Figur 23 Figur 24 
Korn 57.10, D3. Korn 57.10, D 8, mit polarisiertem Licht. 


ZAMP IXa/15 
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getragen. Vergleichshalber sind jeweils auch die zugehérigen Aufnahmen in 
durchscheinendem Licht gezeigt. . 

Die Erfahrung lehrt nun, dass langliche Einkristalle gleich wie Blasenreihen — 
vom Zentrum weg in Richtung des Wachstums zeigen (Figuren 14, 26 und 28). 
Fallt nun eine Diinnschnittebene nicht in die Hauptebene des untersuchten 
Kornes, so ergeben sich schleifende Schnitte mit den Einkristallen, was ein 
Bild gemass Figur 20 ergibt oder, bei noch grésserer Exzentrizitat, zum Teil - 
deutlich zu erkennende Lagen senkrecht zur Wachstumsrichtung (Bildmitte 
der Figur 24). . 


5.2 Einkristallanordnung 1m Wachstumszentrum 


Eine Untersuchung des Aufbaues eines Hagelkornes hat stets von seinem 
Ausgangspunkt, dem sogenannten Wachstumszentrum auszugehen. Vielfach 
ist zwar dessen Lage bereits aus der Anordnung der Luftblasen ersichtlich, die 
Betrachtung des kristallinen Aufbaues erst gibt uns eine zwingende Aussage 
dariiber. 

Figuren 27 und 28 gewahren einen Eindruck iiber das Korn 57.7 (siehe auch 
Figur 1!), dessen Ausgangspartikel eine Frostgraupel kegeligen Aussehens 
bzw. eine konische Reifgraupel war. Das Gefiihl eines Zentrums mit einer 
kugeligen Frostgraupel geben dagegen Figuren 29 und 30, wahrenddem Figuren 
31 und 32 wohl einen urspriinglichen Wassertropfen mit einer nicht ganz ge- 
schmolzenen kegeligen Reifgraupel im Innern zeigen. Ein einkristallartiges 
Wachstumszentrum, basierend auf einer Klareispartikel, ist in den Figuren 33 
und 34 dargestellt. 


5.3 Anordnungen von Einkristallen zweiter und weiterer Generationen 


Die Anordnung der Generationen kann derart sein, dass das Wachstums- 
zentrum in der geometrischen Mitte (Figuren 35 und 36), exzentrisch (Figuren 
37 und 38) oder auch an der Kornoberflache liegen kann (Figuren 39 und 40). 


Die Hauptursache dieser Variation liegt in den speziellen aerodynamischen 
Verhaltnissen beim Wachstum. 


5.4 Einkristallgréssen 


Die Abhangigkeit von der Schnittlage, die Vielfalt der Kornformen und 
Schalen verunméglichen es, die in der Diinnschnittechnik normalerweise zu 
verwendenden Formeln zur Ermittlung der Einkristallgréssen [10] zu gebrau- 
chen. Vergleichbar diirften héchstens mittlere Kornflachen in Hauptschnitten 
sein! Eine Ausnahme hierzu bilden feinkristalline Zonen (zum Beispiel bei 
Korn 57.11, D 2, Figur 38), deren Bild mehr oder weniger unabhangig ist von 
der Exzentrizitaét. Die fiir sie ermittelten Korngréssen diirfen auf Werte 
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Figur 25 Figur 26 
Korn 57.19, D 2. Radial angeordnete Einkristalle, Korn 57.19, 


D 2, mit polarisiertem Licht. 


Figur 27 Figur 28 


Wachstumszentrum, basierend auf konischer Wachstumszentrum, basierend auf konischer Reif- 
Reifgraupel, Korn 57.7, D 2. graupel, Korn 57.7, D 2, mit polarisiertem Licht. 


Figur 30 
Wachstumszentrum, basierend auf kugeliger Frost- 
graupel, Korn 57.10, D 1, mit polarisiertem Licht. 


Figur 29 
Wachstumszentrum, basierend auf kugeliger 
Frostgraupel, Korn 57.10, D 1. 
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Figur 31 Figur 32 
Wachstumszentrum, basierend auf nicht ganz Wachstumszentrum, basierend auf nicht ganz 
geschmolzener Reifgraupel, Korn 57.A, D 1. geschmolzener Reifgraupel, Korn 57.A, D 1, 


mit polarisiertem Licht. 
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Figur 33 Figur 34 
Wachstumszentrum, basierend auf Klareis- Wachstumszentrum, basierend auf Klareispartikel 
partikel Typ II, Korn 57.5, D 2. Typ II, Korn 57.5, D 2, mit polarisiertem Licht. 


Figur 35 Figur 36 
Wachstumszentrum geometrisch zentrisch, 
Korn 57.20, D 1. 


Wachstumszentrum geometrisch zentrisch, 
| Korn 57.20, D 1, mit polarisiertem Licht. 
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zwischen 10-* und 10-1 mm? abgeschatzt werden, was fiir die kleineren den 
Schluss nahe liegen lasst, dass sie durch ebensogrosse Wolkentropfen (~100 uw) 
entstanden sein kénnten. Mitbestimmend hierzu ist auch die Tatsache der vielen 
eingebauten Luftblasen der betreffenden Zonen, die auf eine grosse Gefrier- 
geschwindigkeit ahnlich wie beim Reifgraupelwachstum schliessen lassen. 


5.5 Wirksame Gefrierkerne 


Die oben angefiihrten Hinweise lassen auch die allgemeine Bestimmung der 


_ Zahl der wirksamen Gefrierkerne nur fiir die kleinkristallinen Zonen zu; die 


ermittelte Kerndichte liegt bei 10° Partikeln pro Kubikzentimeter Eis, was Luft- 
gehalten von rund 600 wirksamen Gefrierkernen pro Liter Wolkenluft (Wasser- 
gehalt 5 g/m*) entspricht, und zwar im fortgeschrittensten Stadium der Hagel- 
bildung! 


6. Zusammenhange zwischen Anordnungen von Luftblasen und 
Eiseinkristallen 


6.1 Relationen zwischen Blasen- und Kristailanordnungen 


Vielfach ergibt erst die Synthese der Beobachtungen an Diinnschnitten in 
polarisiertem und durchscheinendem Licht ein ganz zuverlassiges Bild des 
Aufbaues eines Hagelkornes. 

Am deutlichsten kann die Relation anhand der Figuren 41 und 42 demon- 
striert werden, wobei ein Vergleich zeigt, dass dem dichten Klareis pallisaden- 
artige Einkristalle, dem mit vielen Blasen durchsetzten Eis feinkristalline Zo- 
nen entsprechen. Diese Gegentiberstellungen kénnen auch anhand friiher ge- 


. zeigter Figuren gemacht werden. — Solche Aussagen sind aber nur allgemein 


giiltig unter der Voraussetzung, dass geniigend Gefrierkerne und/oder ahnliche 
Gefriergeschwindigkeiten auftreten. Fehlt zum Beispiel die erstgenannte Vor- 
aussetzung, so kann das in Figur 20 gezeigte Bild der zonenweise durch Blasen- 
schalen hindurchdringenden Einkristalle ersichtlich werden. 

Beobachtungen am Korn 57.11 (Figur 44) in polarisiertem Licht wiirden 
auf den ersten Blick keinen Schalenaufbau vermuten lassen, wie er in Figur 43, 
die in durchscheinender Beleuchtung aufgenommen wurde, erkennbar ist. 
Blasenflichen brauchen also nicht unbedingt mit einer feinkristallinen Zone 
gekoppelt zu sein; bei geringer Dicke kénnen sie auch lediglich die Ausgangs- 
flache neuer Pallisadenkristalle bilden. 

Eine Betrachtung der feinkristallinen Zonen in einer etwa 2,5 mm dicken 
Schicht weist auf Ballungen von Lufteinschliissen hin (Figur 45), denen auch 
Gruppierungen von Einkristallen entsprechen (zum Beispiel Figur 46). Da- 
zwischen auftretende Klareiskandle (Figur 47) sind meist nicht als Einkristallite 
gefroren (Figur 48); trotzdem diirfte aber eine értliche Uberlagerung zweier 
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Figur 37 Figur 38 
Wachstumszentrum geometrisch exzentrisch, Wachstumszentrum geometrisch exzentrisch, 
Korn 57abi apie? Korn 57.11, D 2; mit polarisiertem Licht. 


Figur 39 Figur 40 
Wachstumszentrum an der Kornoberflache, Wachstumszentrum an der Kornoberflache, 
Korn 57.2, Korn 57.2, mit polarisiertem Licht. 


Figur 41 


Figur 42 


Korn 57.15, D 1, mit polarisiertem Licht. 


Kor 57:15" Dil: 
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Figur 43 Figur 44 


Konzentrische Kugelschalen Einkristallanordnung im Zentrum des Kornes 
des Kornes 57.11, D1. 57.11, D1, mit polarisiertem Licht. 


es 


Figur 46 


Zusammenballungen von Lufteinschlissen, Einkristallgruppen, Korn 57.B, 
Korn 57.1, D3. mit polarisiertem Licht. 


Figur 47 


Klareiskanale, Korn 57.A, D 1. Entsprechende Einkristallanordnung, 


Korn 57.A, D 1, mit polarisiertem Licht. 
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Wachstumsphasen stattgefunden haben. Es wiederholt sich hier also im Gros- 
sen, was in friiheren Stadien zu Reifgraupeln und nachher zu Frostgraupeln 
fiihrte: Es werden zunachst luftige Schichten aufgebaut, zum Teil in Form von 
Protuberanzen, die spater verdichtet und ausnivelliert werden durch Anlage- 
rung von Wassertropfen bei kleinerer Gefriergeschwindigkeit. 


6.2 Einbauorte der Luftblasen beziiglich der Korngrenzen 


Die Untersuchung, ob der Einbauort der Luftblasen in das Eis abhangig sei 
von den Einkristallen, ergab keine merkbare Bevorzugung oder Benachteili- 
gung der Kornbegrenzungsflachen. 


Fir die Unterstiitzung und Férderung dieser Arbeiten méchte ich vor allem 
Herrn Dr. M. DE QuERVAIN danken, sowie den verschiedenen Hagelkorn- 
sammlern aus der ganzen Schweiz, die durch ihren Einsatz das nétige « Roh- 
material» zur Verfiigung stellten und damit diese Arbeit erst erméglichten. 

Die hier beschriebenen Untersuchungen wurden im Rahmen des For- 
schungsprogrammes der Eidgenéssischen Kommission zum Studium der Hagel- 
bildung und der Hagelabwehr an der Forschungsstelle Weissfluhjoch durchge- 
fiihrt ; zur Finanzierung standen teilweise Gelder des Schweizerischen National- 
fonds zur Verfiigung. 
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Summary 


Photographs and pictures of thin sections in transparent or polarized light 
show a number of different possibilities of-external and internal structure of 
hailstones. The relationship between air bubbles in the ice and the single crystals 


are explained,also in connection with the concentration of the active atmospheric 
ice nuclei. 


(Eingegangen: 21. Januar 1958.) 
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Uber die Darstellung einer Flugbahnschar und ihrer 
St6rungen fiir Feuerleitgerate 


Von Ernst RotH-DEsMEULES, Luzern 


1. Einleitung 


Die vorliegende Untersuchung befasst sich mit einer Darstellung der Flug- 
_bahnen, wie sie in Feuerleitgeraten der terrestrischen Fliegerabwehr Verwen- 
dung finden kann. Aufgabe eines Feuerleitgerdates ist es, aus den gemessenen 
Positionen eines Flugzeuges dessen zukiinftige Bahn zu extrapolieren und ent- 
sprechend die Schiesselemente zu liefern, nach denen die Flabkanone zu richten 
ist. Dies hat zur Voraussetzung, dass ein Rechengeradt vorhanden ist, das die 
_ ganze Flugbahnschar erzeugen kann und ausserdem auch alle Stérungen einer 
Flugbahn zu beriicksichtigen vermag. Da eine laufende Neuberechnung der 
Flugbahn auf Grund der Differentialgleichungen unter Beachtung der momen- 
tanen St6rungen mit den heute zur Verfiigung stehenden Rechenautomaten 
immer noch viel zu zeitraubend ist (abgesehen von deren Umfang und dem 
taktischen Einsatz) und es sich, genau besehen, ausserdem um ein Randwert- 
problem mit komplizierten, nichtlinearen Randbedingungen handelt, das 
einer direkten mathematischen Behandlung noch kaum zuginglich ist, bedeu- 
tet dies, dass mindestens vorlaufig von schon zum voraus berechneten Flug- 
bahnen ausgegangen werden muss. Es muss daher die ganze Flugbahnschar 
im Feuerleitgerat gespeichert vorliegen, damit alle Daten in kiirzester Zeit zur 
Verfiigung stehen, soll der Zweck erfiillt werden kénnen, was bei den heutigen 
_ grossen Flugzeuggeschwindigkeiten bereits ein erhebliches Problem darstellt. 

Im folgenden bleiben alle Fragen, die den Aufbau eines Feuerleitgerates be- 
treffen, insbesondere auch das Problem der Extrapolation des Flugzeugkurses 
unberiihrt. Es wird vorausgesetzt, dass der Treffpunkt genau bekannt sei. Es 
interessiert hier nur die mathematische Form, in der die Flugbahnschar ge- 
speichert werden kann. 


2. Die Flugbahnschar und ihre Parameter 


2.1 Die Flugbahnschar 


Die Flugbahnen seien auf ein rechtwinkliges Koordinatensystem bezogen, 
das in einer vertikalen Ebene (Schussebene) durch das Geschiitz O(0,0) und 
den Treffpunkt T(x, y) liege (vgl. Figur 1). Nun hangt die Flugbahn von einer 
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.grossen Zahl von Parametern ab, zunachst natiirlich von der Zeit t, den An- 
fangsbedingungen vp, é) und weiteren Funktionen, wie der Luftdichte 6, der 
Lufttemperatur 6, dem Wind w1). Die Koordinaten eines Flugbahnpunktes 
hangen damit von den folgenden Gréssen ab: 


4 = x«x(y,b; U9, Eq; 0, 8, w), yes Vi, £205, Bat0,.0,42)., 


T(xy) 


Y. 


Figur 1 


Die Flugbahn im Koordinatensystem. 


Da die Flugbahnen fiir eine Standardatmosphare berechnet werden, ist es 
zweckmassiger, sie durch die Anfangsbedingungen, die St6rungen der Anfangs- 
bedingungen sowie die Abweichungen der wirklichen Atmosphare zu charakte- 


risieren ®) : 


4 = 4, ti Ug, €g; Avg, AE,, 40), y= {V, £3.0_; 29} Lp, Leo, Ole 


In einem Feuerleitgerat sind haufig nicht die rechtwinkligen Koordinaten x, 


des Treffpunktes vorhanden, sondern seine Polarkoordinaten, namlich die Ent- . 


fernung y und der Lagewinkel A. In diesem Falle gilt entsprechend wie oben 
f=r (t; ¥g,A+ 0; Av, Ae; Ad), A=Al(t, 9,4 +0, Avo, Ao, Ad), (2.1.1) 


wo der Schusswinkel o = €) — A ist. 
Die Flugbahnschar hangt unter diesen Voraussetzungen von 5 Parametern 


t, 6, Avy, Ao, Ad ab, da die Anfangsgeschwindigkeit Up als eine Konstante be- 
trachtet werden darf. 


1) Im folgenden soll von den Effekten, die normal zur Flugbahn wirken, wie Derivation und 
Seitenwind, abgesehen werden. Ihre Beriicksichtigung im Feuerleitgerat geschieht ohne Schwierig- 
keiten auf eine entsprechende Weise. 

2) Die Temperaturstérung AO bleibt im weitern unberucksichtigt, obwohl ihr Einfluss nicht 
immer vernachlassigbar ist. Ebenso wird auf den Windeinfluss nicht eingegangen; er kann im 
Prinzip auf eine Stérung der Anfangselemente zuriickgefiihrt werden. 


| 
| 
t 
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2.2 Realisation der Flugbahnschar 


Damit stellt sich die Aufgabe, die Flugbahnschar (2.1.1) in Abhangigkeit 
von ihren Parametern in einer geeigneten Form zu realisieren. Nun ist die 
mechanische oder elektrische Realisation einer Funktion von fiinf Variablen nicht 
ohne weiteres méglich. Handelt es sich um eine Funktion von nur zwei Varia- 


_ blen, etwa ¢ und é9, dann kann sie durch Kurvenkérper erzeugt werden (2, 4] sy. 


Da deren Herstellung nicht nur kostspielig ist, sondern auch ihre Genauigkeit 
beschrankt, so sucht man sie zu vermeiden. 

Bedeutend einfacher ist es, wenn nur Funktionen mit einer einzigen Varia- 
blen auftreten, da diese sowohl mechanisch (Kurvenscheiben, Gelenkwerke [6] 
usw.) als auch elektrisch (Potentiometer [5], Funktionsdrehkondensatoren) 
verhaltnismassig leicht, genau und zugleich robust hergestellt werden kénnen. 
— Dies fiihrt zum Problem, die Funktionen (2.1.1) naherungsweise aus Funk- 
tionen mit einer Variablen zusammenzusetzen, wobei diese Funktionen noch 
durch Addition, Subtraktion und Multiplikation verkniipft werden diirfen. 
Dabei wird man aus Griinden der Betriebssicherheit und der Genauigkeit da- 
nach trachten, mit einem Minimum an Rechenelementen auszukommen. Die 
Aufgabe der Zerlegung in einfache Elemente kénnte etwa durch Versuche ge- 
lést werden, doch stésst dies auf manche Schwierigkeiten, weil die Zusammen- 
hange ziemlich kompliziert sind. Es wurde vorgeschlagen [3], die einzelnen 
Funktionen auf statistischem Weg zu bestimmen. Im folgenden Abschnitt soll 
kurz darauf eingegangen werden. 


2.3 Statistische Approximation. Eine Alternative 


Es soll hier die Grundidee fiir den Fall einer beliebigen*) Funktion F(x, y) 
mit zwei Variablen angegeben werden. Man macht den Ansatz 


F(x, y) = Xf (x) + YE (y) + XP (x) Yo (y) +X3(*) Vey). (2.3.1) 


Die Funktion F(x,y) sei fiir die mn diskreten Werte (x;, yj), (@=1, 2, ..., m; 
j =1, 2,..., 2) gegeben. Man darf dabei voraussetzen, dass der Mittelwertvaller 
gegebenen Funktionswerte null sei (was durch Subtraktion einer Konstanten, 
des Mittelwertes, von den einzelnen Funktionswerten stets erreichbar ist). Die 
Methode der kleinsten Quadrate liefert dann fiir die ersten Glieder 


1 1 ; 
Xf (x) = = Dy Fle, Vi)» YF (y;) = —e D Uiy, 5) 
i] i) 
Sind nun die Reste 
R(G6¢5°V3) =F (43, V3) — [XP (x,) + MA (oy) 


3) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 250. 
4) Es darf sich dabei auch um rein empirisch gegebene Funktionen handeln. 
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dann folgen die nachsten Glieder aus 


LR, V3) °Y; f(y) 
F j 
Xx (%;) = LYFE (ys)? 
7 


und entsprechend Y(y,). Ahnlich kénnen auch die folgenden Glieder X3, YJ _ 
erhalten werden. 

Die Methode lasst sich auf Funktionen mit ae als 2 Variablen erweitern, 
doch wird der Rechenaufwand bald sehr gross, da Summanden mit einer 
Funktion, mit Kombinationen von zwei, drei und mehr Funktionen auftreten. 
Kennt man den Einfluss der verschiedenen Glieder einigermassen, dann asst 
sich der allgemeine Ansatz entsprechend verkiirzen. 

Um zu einem Uberblick iiber den Ansatz zu gelangen, besteht im Falle der 
Flugbahnen die Méglichkeit®), die Tatsache auszuniitzen, dass die Funktionen 
als Lésungen von Differentialgleichungen gegeben sind. Allerdings sind diese 
nichtlinearen Differentialgleichungen 4. Ordnung bei Beriicksichtigung aller 
Einfliisse so kompliziert, dass keine explizite L6sung méglich ist und nur nume- 
rische Verfahren zum Ziel fiithren. Der im folgenden eingeschlagene Weg beruht 
darauf, die Differentialgleichungen durch approximative zu ersetzen, die eine 
Lésung durch bekannte Funktionen gestatten, so dass deren Bau leicht zu — 
tiberblicken ist und Schliisse tiber die Form des Ansatzes (2.3.1) gezogen wer- 
den kénnen®). Dabei ergibt sich sofort noch eine Zerlegung des Problems, da die 
Normalflugbahnen und die Stérungen getrennt behandelt werden diirfen, 
letztere mit Hilfe der St6érungstheorie. — Die Erfahrung zeigt, dass man auf 
diesem Wege zu Hinweisen gelangt, die den noch erforderlichen Rechenauf- 
wand erheblich verringern. 


2.4 Die Differentialgleichungen der Bewegung 


Werden nur der Luftwiderstand und das Geschossgewicht beriicksichtigt, 
dann lauten die Differentialgleichungen der ebenen Geschossbewegung [1] 


f= WL yee Ws Se (2.4.1) 


v / Vv oe 


wo v?= x+y? und W die Widerstandsfunktion ist. Vernachlassigt man die 
Abhangigkeit von der Machschen Zahl und wahlt fiir die Luftdichte einen ex- 


°) Die andere Méglichkeit, etwa Potenzreihen oder Fourier-Reihen zu Hilfe zu nehmen, 
scheitert gewohnlich an der schlechten Konvergenz fiir grosse Zeiten f. 


6 ee : 
) Dieser Weg ist offenbar immer dann gangbar, wenn die zu realisierenden Funktionen 
durch Differentialgleichungen bestimmt sind. 
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ponentiellen Verlauf, dann gilt fiir W der Ansatz 
Wee 0 0ri gy", (2.4.2) 


wo c, a und” Konstanten bedeuten. 
Fiir die Anwendungen bei Feuerleitgeraten ware es gegeben, die Differen- 
_tialgleichungen in Polarkoordinaten zu formulieren, um direkt die gewiinschten 
Gréssen zu erhalten. Die entsprechenden Gleichungen sind aber wenig itiber- 
sichtlich. 


3. Darstellung der Normalflugbahn 


3.1 Lésung der Differentialgleichungen 


Die Schwierigkeiten bei der expliziten Lésung des Systems (2.4.1) legen 
in der Gestalt der Luftwiderstandsfunktion (2.4.2). Fiir a + 0 fiithrt allein der 
Fall = 1 auf eine lineare Differentialgleichung, deren Losungen aber insofern 
nicht explizit angebbar sind, als Quadraturen stehenbleiben. Fiir a = 0 ist nur 
die Geschwindigkeit explizit ausdriickbar, wahrend die Koordinaten wieder 
durch Quadraturen zu bestimmen sind. Einzig der Fall n =1 fihrt ans Ziel 
und soll daher im folgenden allein betrachtet werden’). 

Durch elementare Rechnungen findet man fiir a = 0 und m =1 unter Be- 
riicksichtigung der Anfangsbedingungen (¢ = 0) 


g=y =O, i == "Uy COS £4 y = Uo SiN Ey 
die Lésung 
1 C—tge 
eer CS on 
mit 
C= item tge 
C Ug COS Eg e 


und daraus fiir die Koordinaten 


Up COS E on, 8 Uo SiN Eq g a ee 
xm LOCOS Po (4 — ¢-et), y=| : +S] 0 ent) — Eg, (3.1.1) 


3.2 Ansatz fiir den Schusswinkel a(t, A) 


Entsprechend den Uberlegungen in Abschnitt 2.3 soll der Schusswinkel 
a(t, A) durch Funktionen 7,°(¢), A7(A) (7 =1, 2, ...) einer eizigen Variablen ¢ 


7) Diese Vereinfachung ist insofern gerechtfertigt, als im Falle von Uberschallgeschwindig- 
keiten der Ansatz W=cv! fiir die gestreckten Flabbahnen eine sehr gute Naherung darstellt. 
Somit ist zu erwarten, dass sich auch iiber den Aufbau des Ansatzes (2.3.1) weitgehende Aussagen 
machen lassen, was durch die Erfahrung bestatigt wird. 
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bzw. A in der Form 
a(t, A) = T7() + AA) + TY - Aga) + --- 


dargestellt werden, wobei natiirlich so wenig Glieder als méglich zu verwenden ~ 
sind. 
Aus den Beziehungen (3.1.1) folgt fiir den Lagewinkel A . 


t 1 es 
tg = tee, see & - Dy COSE, 


Beriicksichtigt man ¢)=A-+ 0 und beachtet, dass o im allgemeinen klein ist, 
dann ergibt sich die fiir nicht zu grosse Zeiten giiltige Naherung 


o=Ty(i) cosey mit Ty) =—%-(1-;4 =). 6.21) | 


4 CU 1 —e—¢ 
In manchen Fallen wird man cos é) + os A setzen diirfen§) : 
o = TY (t) cosa. (3:22) 
Geht man in der Naherung einen Schritt weiter, dann ist 


o [1+ T(t) sind] = T(t) cosa 
oder auch 
o = Ty’(t) cosA — TY (t) sin2 A » - (3.273) 


mit T(t) = [Ty (0)}?/2. 
Die Darstellung (3.2.3) zeigt, dass die ersten Glieder Ty’, A? gar nicht er- 
forderlich sind®). 


3.3 Ansatz fiir die Entfernung r(t, A) 


Fir die Entfernung 7 des Treffpunktes soll in analoger Weise gelten 


nt, A) = Ti (0) + AL(A) + TI) - AEA) + 


8) Diese Naherung schliesst man bereits aus der Vakyumbahn, ftir die 


2 


29 l : : 
r= > — — cos(A +o) sino 
g cos? A 
gilt, woraus 
ry 
ja 208 cosdA = T§(t) cos 


folgt. 
9) In der Folge werden die trigonometrischen Funktionen stehengelassen, obwohl es haufig 
zweckmassig sein wird, sie durch leicht davon abweichende Funktionen A,(A) zu ersetzen. 


; 
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_ 


y kann aus x oder y berechnet werden. Man unterscheidet zweckmassig zwei 
' Faille 


1.* >y (untere Winkelgruppe); 2. x < y (obere Winkelgruppe). 


Im ersten Falle eignet sich fiir 7 die Darstellung 


x 


U 
= seep =o (l—e-%) cosa — 2 (1—e-*") sino tgd. (3.3.1) 


Bedenkt man die Gréssenordnung von o sowie die Approximation (3.2.2), so 
erhalt man 


jim “9 (1 —eret) — (1 — e-¢?) Tye(é) sind, 


woraus man die Form 
r= T17(t) — TZ) sind (SP5ea) 


abliest, mit den Abkiirzungen 
RWS l-e, LOO 120. 


Um die Naherung einen Schritt weiterzutreiben, kénnen in (3.3.1) die ersten 
Glieder der Reihenentwicklungen von coso und sino eingefiihrt werden. Dies 
liefert 

r(t, A) =T,*" (2) —Tj(é sina — Tj (t) cos2 A, (373.3) 


_ wo sich T,*’ und Tj durch Tj und 7,’ ausdriicken lassen, womit gegebenenfalls 
‘noch Funktionselemente eingespart werden kénnen?°), Jedenfalls lehrt die 
Darstellung (3.3.3) sofort, dass eine Funktion Aj(A) nicht auftritt. 
Fiir den Fall des obern Winkelbereiches ergibt sich 7 aus 7 = y/sin A. In einer 
ersten Naherung folgt wieder die Beziehung (3.3.2) mit denselben Funktionen. 
Entsprechend kann man die Darstellung um ein weiteres Glied Tj -Aj er- 


ganzen. 


4. Ansatz fiir die Stérungen der Flugbahn 


4.1 Die allgemeinen Stérungsdifferentialgleichungen 


Die Flugbahnen sind fiir Normalbedingungen berechnet, von denen die 
wirklichen Verhiltnisse mehr oder weniger stark abweichen kénnen. Solche 
Stérungen sind Anderungen des Elevationswinkels ¢,, der Anfangsgeschwin- 
digkeit v), Abweichungen der Luftdichte vom Standardverlauf, Wind usw. 


10) Sollte noch ein Ausdruck mit cos A erforderlich sein, dann kann dieser mit T3sinA zu 
einem Glied 73” sin(A — A») vereinigt werden. 
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Sind die Stérungen klein, dann kénnen die Abweichungen des Ortes vom 
Normalort in erster Naherung mit Hilfe der Stérungstheorie berechnet werden. 

Es seien zundchst nur die Stérungen von c und v, betrachtet, und zwar 
seien die neuen Werte c + y und v, + Avy. Fiir die gestérten Koordinaten der 
Flugbahnpunkte gelte 


X=xz+6&, Y=y+y%, (4.1.1) 


wo & und 7 die Verschiebungen des Flugbahnpunktes bedeuten. Nun erfiillen 
X und Y auch die Differentialgleichungen (2.4.1). Mit Hilfe des Luftwider- 
standsgesetzes (2.4.2) kénnen beispielsweise durch Reihenentwicklung die zu- 
gehorigen Variationsgleichungen abgeleitet werden. Nach einer langeren 
Rechnung folgen fiir € und 7 die linearen Differentialgleichungen 


ga ceo 1 ti] e+ 42 hadnt 2a, 

UV U c 

ae ° (4.1.2) 
j= ee orn 42 E+ [m0 41] pans 25h 


Diese Differentialgleichungen’ in die x und y als Lésungen von (2.4.1) eingehen, 
lassen keine explizite Integration durch bekannte Funktionen zu, so dass man 
auf numerische Methoden angewiesen ist. 

Selbst im Falle n = 1, fiir den sich die Gleichungen erheblich vereinfachen, 
ist diese Vereinfachung nicht wesentlich. Aus diesem Grunde soll weiterhin nur 
der Spezialfall 7 =1, a = 0 betrachtet werden, um explizite Lésungen zu er- 
halten. Auf diese Weise kann die Form der zu approximierenden Funktionen 
genauer untersucht werden. 


4.2 Explizite Losung fiir den Fall n=1, a=0 


Obwohl in diesem Falle fiir die Normalflugbahn eine konstante Luftdichte 


6 vorausgesetzt wird, lassen sich Dichtestérungen A6 beriicksichtigen, dazufolge 
der Beziehung 


die Dichtestérungen als Stérungen des ballistischen Koeffizienten c aufgefasst 
werden kénnen. Es sollen die folgenden. Stérungen berechnet werden: 

1, Stérung der Anfangsgeschwindigkeit v)—> vy + Av9; 

2. Stérung des Elevationswinkels ¢)—> @) + Ae; 


3. Stérung des ballistischen Koeffizienten c>c +. 
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Die Variationsgleichungen (4.1.2) nehmen nun die einfache Gestalt an: 
Etc&=—y#, hteoq=—yy. (4.2.1) 


Die inhomogenen Glieder mit % und y sind nach (3.1.1) bekannte Funktionen 
der Zeit. Die drei erwahnten Falle unterscheiden sich durch die Anfangsbe- 
dingungen und die inhomogenen Glieder. 


4.2.1 v9-Stérung. Das System (4.2.1) ist mit y = 0 unter Beriicksichtigung 
der Anfangsbedingungen fiir ¢ = 0 


E(0) = 7(0) = 0, E(0) = Av, COS Ep, 7)(0) = Ady sine, 


zu lésen. Es ergibt sich 


COS Ey 


See ta Avy | 6.2.2) 


f= (1 =~) Av, Ng = 


4.2.2 &5-Storung. In-diesem Falle ist y=0, und die Anfangsbedingungen 
lauten (¢ = 0) 
£(0) =7(0) =0, —€(0) =—upsiney dey, (0) = + 19 cose4 Aey. 


Die Lésung ist 


Up SIN Eg ( 
c 


Ba eee eagle (4.9.3) 


—é = — 1—e°) Ae, Ne, = + 


4.2.3 c-Stérung. Die Differentialgleichungen (4.2.1) lauten, ausfiihrlich 
geschrieben, 


zt : Bp 5 i 2 ses g 
Etc&=—yucose_e”, nten=-y (2 + Up sin e9) e Cae 
Anfangsbedingungen fiir ¢ = 0 
£0) (0) =0, £0) = (0) = 0. 


Als Lésung findet man 


§,= — SS [et ete} y, 
as eee fie ote (4.2.4) 


1 & fore? —20—e"MN 


ZAMP IXa/16 
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4.2.4 Folgerungen. Aus den Lésungen (4.2.2) und (4.2.3) ergeben sich ohne ; 
weiteres zwei bemerkenswerte Beziehungen fiir die durch die Stérungen Av) 
und Ae, bewirkten Verschiebungen des Treffpunktes T. 


Satz 1. Ist der Luftwiderstand proportional zur Geschwindigkeit und die 
Luftdichte konstant, dann ist die durch eine Stérung Avy der Anfangsgeschwindig- 
keit bewirkte Verschiebung des Treffpunktes parallel zur Abschussrichtung. 


Satz 2. Unter denselben Voraussetzungen wie in Satz 1 ist die durch eine 
Storung Ae, des Elevationswinkels bedingte V erschiebung des Treffpunktes normal 
zur Abschussvorrichtung. 


Naherungsweise gelten diese Eigenschaften der Treffpunktverschiebung 
auch fiir die allgemeinen (empirischen) Luftwiderstandsgesetze!). 


4.3 Die Korrekturen im Falle einer vy-Storung 


4.3.1 Die Korrekturen At,, und Ae,,. Unter dem Einfluss einer Stérung 
gelangt das Geschoss in der Zeit ¢ nicht nach dem Punkt T(x, y), sondern nach 


TO , Y). Soll es wiederum durch den Punkt T gehen, dann muss es unter dem 
Winkel ¢) + Ag) abgefeuert werden und wird dabei die Zeit ¢ + At bendtigen 


Yi 


ie 


Figur 2 
Die gestorte und die korrigierte Flugbahn. 


(vgl. Figur 2). Zu einer gegebenen Stérung Av, sind nun die Korrekturen Ag, 
und A?,, zu berechnen. Es gilt in erster Naherung 


—f=XAt, —y=YAt, (4.3.1) 


1) Man uberzeugt sich leicht, dass die Satze 1 und 2 auch fiir die Flugbahn im Vakuum gelten. 
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_ wobei die Stérungen & und 7 durch die folgenden Ausdriicke gegeben sind: 
Crees # ape ie. (4.3.2) 
und weiter ist . 
X=44+8,46,, Y=pt, + te: (4.3.3) 
Werden die Beziehungen (4. 2. 2) ‘und (4.2.3) unter Beachtung von (4.3.2) und 


(4.3.3) in die Gleichungen (4.3.1) eingesetzt, dann bestimmen sich die unbe- 
kannten Gréssen Af, und Ae, aus dem folgenden Gleichungssystem: 


oS EO (1 — 07!) Aey + (Up + Avy) cosey e~“ At 
— vp sine, e*' Ae, At = ae (1—e~“) Avy, 
| (4.3.4) 
= Hoe. (1 —e~°") Aeg+ [(o + Avg) sine, ec — = i= e-*)] At 
+ J COS&, € *' As, Ai = pas (1 — e=) Avy. 


Da Ae, und Af kleine Gréssen sind, diirfen — wie eine Abschatzung zeigt — in 
(4.3.4) die quadratischen Glieder ohne nennenswerten Fehler weggelassen wer- 
Poden2). 
Als Lésung findet man die folgenden Korrekturen 


= (l-en" 
At, =— z Avy (4.3.5) 


% g SiN Ey 
. a ct 
UgiGmiae : (l—e-*%) 


und 
& COS £4 = (Leet) 
Lg z Avy. (4.3.6) 


% g SiN €y 
= —ct 
Dee oh oes (1—e-¢?) 


Man erkennt sofort, dass zwischen den beiden Werten die einfache Relation 


ewe: a COS Ey At, (4.3.7) 


0 


besteht. Sie kann dazu beniitzt werden, um aus der Kenntnis der einen Kor- 
rektur sofort auf einfache Weise die andere Korrektur zu berechnen. 


12) Dies lauft darauf hinaus, im System (4.3.1) X durch # und ebenso Y durch y zu ersetzen, 
was im allgemeinen eine zulassige Approximation ist. 
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4.3.2 Approximation durch Funktionen mit einer Variablen. Um nun die 
Korrekturen analog wie die Flugbahnelemente durch Entwicklungen von der ; 
Form (2.3.1) darzustellen, werde (4.3.5) in eine Reihe entwickelt. Es wird 


I te Po (1 — en“) < [1 4 ERB (et —1) + ---| 
0 


CU, ahs 


Unter Verwendung der Abkiirzungen 


ct 
T;**(t) = g 1 T;*(t) = Bd da (e2** — e°') 


Oy vec 
kann man dafiir 
Av, 


Vo 


At, = —Ty (t) (T(t) + T7*(6) siney + ++] 


‘ (4.3.8) 
setzen. Fiihrt man anstelle des Elevationswinkels ¢, den Lagewinkel 4 ein, so 
ergibt sich die weitere Darstellung 


At, = —Ty(t) (LF"(@) + Til) sind + T*%(6) cos24] 22 (4.3.9) 
M 0 


mit den Funktionen 


1 o U, il vu, o 
TA*() =T) +> RO): LO, Tt) = 5 I) - TO. 
In manchen Fallen wird die einfachere Form 


At, = —Ty() (T°) + T"@) sina] (4.3.8*) 
0 
ausreichen. 
Wie die Darstellungen (4.3.8) und (4.3.9) zeigen, lasst sich die friiher in 
(3.3.2) schon aufgetretene Funktion T,(¢) als Faktor abspalten, was natiirlich 
eine Vereinfachung bedeutet. 


Die Korrektur des Schusswinkels o folgt aus (4.3.8) sofort mit Hilfe der 
Beziehung (4.3.7) 


Ao, = —Ty(t) a COS €y [7;°(t) + Ty(t) sineg ++: oe : 


(4.3.10) 


4.3.3 Der Fall grosser Stérungen Avg. Da in der Praxis oft Stérungen auftre- 
ten, fiir die die lineare Stérungstheorie nicht mehr ausreichend ist, so konnen in 


einfacher Art und im allgemeinen geniigend genau noch in Av, quadratische 
Glieder hinzugefiigt werden, indem man setzt 


At,, = — Ty (T(t) +++] 6 (2) 


U9 
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mit 


© (<2) SL dl +a (<2), (4.3.11) 


U9 Vo Gy 
wo « eine Konstante bedeutet. 


4.3.4 Die Distanzkorrektur Ar, . Unter Umstanden ist es vorteilhaft, auch 
die Korrektur der Zielentfernung einzufiihren 


y* (t+At, A) =r(t, 4) + Ar, (t, a). 


Hier bedeutet 7* die gemessene Distanz. Daraus lasst sich direkt die korrigierte 


Flugzeit entnehmen. 
Um Ar,, zu erhalten, geht man am besten von der Beziehung 


A65 tee BS 


ed ~~ 
UY cos/ COS Eg 


aus, wo wegen (4.2.2) und (4.2.3) 


Av, 


E=F, + & =coséy Ty (it) - ars sin &9 Ty (t) Ago 
0 


zu setzen ist. Mit (4.3.10) folgt schliesslich 


Ar, = Tit) 1+ Ti") + Fe" sind Ae, (4.3.12) 
0 
Es bedeutet 
TAj=£ TOT) umd 7") = OT). 


Entwickelt man in den Ausdriicken auf der rechten Seite noch die Exponential- 
funktionen nach der Zeit, so erkennt man, dass die Reihen erst mit /+ beginnen, 
so dass meist mit geniigender Genauigkeit 


F Av 
Ar, = Frit) [1+ = t+ J, t sine] A“ (4.3.13) 


Uo 0 


gesetzt werden darf. 
Auch hier lasst sich natiirlich fiir den Fall grosser Geschwindigkeitsande- 


rungen Av, als Faktor auf der rechten Seite ein quadratischer Ausdruck 
G@=—*+a ( 


setzen. Erwahnenswert ist, dass in den obigen Darstellungen wiederum der 
Faktor 7/(é) auftritt. 
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4.4 Die Korrekturen im Falle einer c-Stérung 


4.4.1 Die Korrekturen At, und Ae,. Fiir die zu einer gegebenen Anderung y 


des ballistischen Koeffizienten c erforderlichen Korrekturen At, und Ae, gel- — 


ten die zu (4.3.1) analogen Beziehungen 


Za Slate Abell aap ae Volts (4.4.1) 


wobei in diesem Falle 


(4.4.2) 


bab ee, N=, +N,» | 
X=¢4+E,46, Yaptaitn, 


zu setzen ist. 

Fiihrt man die Ausdriicke (4.2.3) und (4.2.4) ein, dann ergibt sich ein 
Gleichungssystem, das. nach der Unterdriickung der kleinen quadratischen 
Glieder die Form 


A sin &9 Ae, — a Cos €q At, = Bcosey 
A cos éq Ae, + (2 SiN Eg — = A) At, = —(B sineg+ G) y | 
0 


annimmt, mit Koeffizienten A, B, G und a, die Funktionen der Zeit sind. Nach 
einer langeren Rechnung findet man fiir die Zeitkorrektur 


By (loti seein Sere [et (e-e# +1) —2 (1—e-)] 
At, ar g ae y (4. 4, 4) 
enact eF 0 (1—e-°t) : 


und fiir die Korrektur des Schusswinkels 


_ § COS Ey et se t— act 
ae C7 Ug ar g sing, r* (4.4.5) 
Sere) |S (e*t—1)] 
Uy C 


Fiir den letzten Ausdruck lasst sich auch 


Mates 2 g COSEy sinhct—ct 
y c2u, | sinhcé#—(1+2h) (coshe#—1) “ 


schreiben, wo k = g sinég/v) c bedeutet. 


4.4.2 Approximation durch Funktionen einer Variablen. Die Entwicklung 
des Ausdruckes (4.4.5) nach Potenzen von ¢ liefert 


rea (oe 
= #4 2 (15 H#— 7) #4 | y (4.4.6) 


on fe 
mt te) 2 
Ae, = COS Eg U lL; ¢2 + 


mit T=14+ 2%, 


| (4.4.3) 
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Bemerkenswert ist, dass die Reihe erst mit #2 beginnt. Konkrete Beispiele 
_ lehren, dass bereits das erste Glied eine sehr gute Naherung darstellt, so dass 
_ man in den meisten Fallen fiir die Schusswinkelkorrektur 
1 Ae, = hy t? cose, © (4.4.7) 
setzen darf. 
_ 4.4.3 Die Distanzkorrektur Ar,. Die Distanzkorrektur ergibt sich aus 


fa gh wie eM 


(2 .GOsd COS Ey © 


Setzt man (4.4.7) in (4.2.3) ein, dann folgt zusammen mit der Beziehung 
(4.2.4) die Formel 


Ar, = —Ty (0) (T7(@) + hy # sine] +. (4.4.8) 


mers ist 7’ =ct/(e* —1)*}-und T, = kp t?. 
Die Formel (4.4.8) zeigt auch diesmal wieder das Auftreten des Faktors 
1 (0). 
* Sind die Stérungen y/c gross, wie dies gelegentlich auftreten kann, dann 
_ lasst sich dies durch die Einfiihrung eines quadratischen Ausdruckes 


paler Acs 


c 


beriicksichtigen; damit wird eine Naherung erzielt, die praktisch ausreichend 
ist. 
5. Ergebnis 


Beniitzt man die Differentialgleichungen der Geschossbewegung mit einem 
naherungsweise giiltigen Luftwiderstandsgesetz, das eine explizite Lésung ge- 
stattet, dann lassen sich daraus wichtige Hinweise gewinnen fiir die zum Auf- 
bau von Flugbahnrechnern in Feuerleitgeraten bendtigten Darstellungen der 
verschiedenen Gréssen durch Funktionen einer einzigen Variablen. Fiir die 
wichtigsten Bahnelemente ergeben sich dabei die folgenden Formeln, die auch 
die Fille grosser, nicht mehr linearer Stérungen erfassen °) 


r(t, A) = Tj (t) — Tj (t) sina — Ty (t) cos2 A, 


*} 


ri) 


Ay, = Fit) 1+ T2"() + TY) sind] d( 


Uo 


Ar, = —Ty(t) [Ti (t) + he #2 sin] ¥ (2), 


; ct 
18) Dieser Ausdruck bleibt auch fiir > 0 endlich, da ran eet =1. 


250 Ernst RorH-DESMEULES ZAMP 


o(t, A) = T,f(t) cosA — T,7() sin2 A, 


: A 
Ao,, = — Tj (t) cose [Ti"(t) + T(t) sin eo] o( =) ; 


Uy %o 
Die einzelnen Funktionen Tj, A’, ..., die durch spezielle Funktionselemente 


zu realisieren sind, kénnen nun auf Grund von exakten Flugbahnberechnungen 
noch naher bestimmt und den konkreten Verhaltnissen angepasst werden. Auf 
diese Weise wird eine Genauigkeit der Darstellung erzielt, die den Anforderun- 
gen der Praxis gerecht zu werden vermag. 

Die vorliegende Arbeit stellt die Erweiterung einer Untersuchung dar, die 
seinerzeit fiir die Firma Contraves AG, Ziirich, ausgefiihrt wurde. Es sei auch 
an dieser Stelle Herrn Dr. A. GERBER, Direktor der Werkzeugmaschinenfabrik 
Oerlikon, Biihrle & Co., und der Contraves AG, Ziirich, fiir sein Interesse an 
dieser Untersuchung bestens gedankt. 
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Summary 


It is shown how the trajectories of antiaircraft shells and their perturbations 
can be approximated with sufficient accuracy in terms containing functions of 
only one single variable. In directors, these functions can be realized by suitable 
function producing calculation elements. 


(Eingegangen: 26. November 1957.) 
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Méglichkeiten und Grenzen des Infrarotsehens 


Von ToNI SCHNEIDER und PETER ScHRANZ, Ziirich 1) 


1. Einleitung 


Bildwandlergerate vollziehen die Transformation von einem Gegenstand in 
ein Bild mit anderem Spektrum und anderer Intensitat durch einen Umwand- 
lungs- sowie Verstaérkungsvorgang in einer unmerklichen Zeitspanne. (Unter 
Gegenstand und Bild wird dabei eine ebene Konfiguration von Strahlungs- 
dichten verstanden.) 

Die Bedeutung solcher Gerate liegt in folgendem: 

a) Man kann Bildwandlerphotokathoden herstellen, die auch auf unsichtbare 

Bereiche des Spektrums ansprechen; 

b) man kann schwache «Gegenstande» auf ein Helligkeitsniveau heben, das 
den von der Leuchtdichte abhangigen Sehfunktionen besser angepasst ist. 


Ausgehend von den Strahlungsgesetzen bei optischer Abbildung, dem Um- 
wandlungs- und Verstaérkungsvorgang in der Bildwandlerréhre sowie den 
Strahlungsdichteschwankungen, behandelt dieser Bericht Grenzen und Még- 
lichkeiten des Infrarotsehens unter der Voraussetzung einer bestimmten Seh- 
anordnung. 


2. Begriffe und schematische Anordnung 


Aus der Fiille der Méglichkeiten wird den folgenden Betrachtungen eine 
Geratevariante zugrunde gelegt, deren schematische Anordnung aus der Figur 
hervorgeht. 


inet, 
‘> D 
BSD, 


<} 
q 


——_— ff —— 


as ta 3 > tt 
Objektiv  Photokathode —_Leuchtschirm Lupe Auge 


1) Physikalisches Laboratorium der Albiswerk Ziirich AG, 


2a2 
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Gesichtsfeld des Infrarotsichtgerates (gesamte auf die wirksame Photo- 
kathode abbildbare Flache) ; 

gesamte wirksame Flache der Photokathode; 

Quadratchen im Gesichtsfeld ; 

optisches Kathodenbild von h?; 

Leuchtschirmbild von h’?; 

Strahlungsdichte von h? [W/cm?]; 

Bestrahlungsstarke’ des optischen Bildes h’? [W/cm?] ; 
Strahlungsdichte des Leuchtschirmbildes h”? [W/cm?]; 
Eintrittspupille des Objektivs [cm]; 

Objektivbrennweite [cm]; 

Transmission des Objektivs [%]; 

effektive Offnung des Objektivs; 1/K; = VT Dol fos 

Abstand Gegenstand—Objektiv ; 

AbbildungsmaBstab der Réhre; 

wirksamer Kathodendurchmesser [cm]; 

Unscharfe des Gerates, bezogen auf die Kathode [cm]; 

relative spektrale Empfindlichkeit der Photokathode [°%]; 
relatives spektrales Emissionsvermégen des Leuchtschirms [°%]; 
relative, spektrale Hellempfindlichkeit des Auges [%]; 

integrale Kathodenempfindlichkeit {[A/W]; 

Schirmausbeute [%]; 

Konstante; hangt von der Schirmstruktur ab; 

AU=U—Us,; 

Beschleunigungsspannung [V] ; 

«dead voltage» [V]; 

Eintrittspupille des Auges [cm]; Funktion der Schirmleuchtdichte [1]?) ; 
Lupenbrennweite [cm]; 

Transmission der Lupe [%]; 

Lupenvergrosserung 


lr, = = : (2.1) 
Gesamtvergrésserung des Geriates 
ee 
ws My 5 (Boe) 


Winkel, unter dem /?2 erscheint, 


(«) = [min]; (2.3) 


eae ee : ‘ : ; A 
) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 259. 
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%¢  Gesichtsfeldwinkel des Gerites 
XK ra Dy, 


a? J? aus geometrisch-optischen Beziehungen folgt: 


P2=L2 2 (565750 = Ve (2.5) 


3. «Sehvorgang¢ » 


Bekanntlich besteht bei optischer Abbildung zwischen der Strahlungs- 
dichte des «Gegenstandes» und der Bestrahlungsstarke des «Bildes» die Be- 
ziehung 


(3.1) 


_wenn nur der AbbildungsmaBstab m< 1 ist. Die Bestrahlungsstarke auf der 
Photokathode, hervorgerufen von einem «Gegenstand», betragt demnach 


nS 


B= } Dee 
oi 8.2) 


unter S wird dabei der von der Photokathode verwertbare Anteil verstanden, 
das heisst 


Si if S(a) KU) da. (3.3) 
S(A) spektrale Emission des Gegenstandes [W/cm*] . 


Die Strahlungsdichte des Leuchtschirmes, hervorgerufen durch E”, ist: 
nS" = E' m,? n = E' mR” ex AU? a,. (3.4) 
7 Umwandlungsfaktor [W/W] . 
Offensichtlich betragt die Strahlungsdichteverstarkung der mit einem 
Objektiv bewaffneten Bildwandlerréhre 


Sie n if ex AUPa, 1 


Bo: Stee Ae a5 
S ~ me 4 Kt me, 4K%. a) 


Im menschlichen Auge gelangt von der Strahlungsdichte S” des Leuchtschir- 


mes nur der Bruchteil 
B= f S” [sb] (3.6) 
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zur Verwertung, wobei 
1 f- Wa 3 
Bo fimvada MM 
M Mechanisches Lichtaquivalent {lm/W] . 


Auf Grund der Dualnatur des Lichtes kann die Abstrahlung eines Leucht- 
schirmflachenelementes als Photonenemission aufgefasst werden. Die ausge- 
sandten Photonen sind iiber die Zeit in solcher Weise verteilt, dass ihre durch- 
schnittliche Zahl pro Flachen- und Zeiteinheit in einem gegebenen Raumwinkel 
proportional der Leuchtdichte ist: 


N=xBh". (3.8) 


N Anzahl Photonen pro Sekunde; x in den Halbraum ausgesandte Photonen- 
zahl pro Lumen. 


Von diesen N in den Halbraum emittierten Photonen kann das mit einer 
Lupe bewaffnete Auge nur den Bruchteil 
2 


auffangen. 

Das fragliche Leuchtschirmelement (h”?) zu sehen, besteht in der Fahigkeit 
des Auges, aus der Menge und Verteilung der auf die Netzhaut auftreffenden 
Photonen Informationen zu schépfen. 

Sei h"? ein h"® benachbartes Leuchtschirmelement, hervorgerufen durch 
zwei benachbarte Flachenelemente in der Gegenstandsebene, dann betragt die 
ihnen zugeordnete Photonenemission 


Not=2#B,ho?t bzw. Nt=xBh"?t. (3.10) 


Da sich das Auftreffen der Photonen aus zufalligen Ereignissen zusammensetzt, 
ist dem Mittelwert der wahrend einer gewissen Zeit die Netzhaut treffenden 
Photonen eine Schwankung tiberlagert, deren quadratischer Mittelwert propor- 
tional zur Quadratwurzel der auftreffenden Photonen ist. Wenn nun die Dif- 
ferenz der von der Netzhaut verwertbaren Photonen nicht grésser ist als die 
Schwankung um N, so ist es prinzipiell unméglich, die fraglichen Flachen- 
elemente (i9*, h*) als getrennt zu empfinden; eine Trennempfindung ist dann 
gewahrleistet, wenn 

N’Ot—-NiOt=KyN'O: (3.11) 
ist, wobei 


@ Quantenausbeute des Auges; Funktion der Schirmleuchtdichte [2]. 


Vol. [Xa, 1958 Moéglichkeiten und Grenzen des Infrarotsehens 255 


_t Speicherzeit; sie ist die Zeit, in der die Photonen in einen Helligkeitsein- 
druck umgesetzt werden. Die Tragheit von Bildwandlergeraten ist aber 
im allgemeinen nicht durch das Auge, sondern den Leuchtschirm gegeben; 

_ das heisst, ¢ muss als Speicherzeit des Leuchtschirmes interpretiert werden. 

& gibt an, wievielmal sich die absorbierte Quantendifferenz aus dem statisti- 
schen Rauschen herausheben muss, damit eine Trennempfindung der beiden 
Quadratchen zustande kommt. Analog oben (siehe unter Speicherzeit) muss 
das fiir die ganze Sehvorrichtung spezifische K gesetzt werden. 


In der Beziehung (3.11) wurde die Untergrundaufhellung des Schirmes noch 
nicht beriicksichtigt ; sie kommt im wesentlichen durch Superposition folgender 
Effekte zustande: 

a) Durch Streustrahlung aus dem Sichtraume; 

b) durch Reflexion in den optischen Teilen ; 

c) durch thermische Elektronenemission der Photokathode sowie der mit 
einem Alkalimetallfilm belegten Oberflachen des Réhreninnern ; 

d) durch Strahlung, die infolge der Transparenz der Photokathode ins Rohren- 
innere gelangt, dort reflektiert wird und an der Kathode Elektronen ablést. 


Offenbar setzt sich die Untergrundaufhellung aus je einem durch das Sicht- 
gerat sowie den Sichtraum verursachten Anteil zusammen. 


7] 
By + Bo=B (amex * Sot Sc). (3.12) 


S, Strahlungsdichte der Streustrahlung aus dem Sichtraum; S¢ gerateseitig 


bedingte Aufhellung. 
Somit erganzt sich (3.11) zu 


N’'Ot—-NOt=KYN'Ot+NiOt+N,Ot. (3.13) 


Dadurch, dass jeder Photonendifferenz eine Leuchtdichtedifferenz und damit 
ein Kontrast entspricht, kann (3.13) in eine Grenzbeziehung umgeformt wer- 
den, die besagt, dass bei gegebenen Geratedaten zu jeder Strahlungsdichte des 
Gegenstandes eine Wertegruppe von Schwellenkontrast, Unschiarfe und Spei- 
cherzeit gehért, die nicht unterschritten werden kann. 


Ee GLI NS ny ire nt Nat NG 
( ae ) Cee Sa (3.14) 
Mit 
2 k2 
, ; ; x Jag = 2 n 3 caf ‘ee Mae T 
fee ut Nc 4G (seas) ore aya Sh he (s AL a eae 
und 


annie) 26.» \P 2 y| aS 
eer Geo-a0} a rs ag Me 
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aus (2.2), (2.3), (2.5), (3.4), (3.9), (3.12) kann (3.14) auch in der Form 


T 1 
2 L n 2 y2 ¢ 
Bos 42 m?, mr K7 site 
Se ASA 
v 
ite é (3.15) 
4 (360-60\? 1 
> kK2 
aes a 2a ) D4, 0 | 
oder 
Fo Nae ak te ali ek 
~ Chat = Goncs ala 
nt B 1 (3.16) 
* ni, K? m2, K%, 
Ss 
dike 3 


geschrieben werden. J’, ein Mass fiir die Grésse (x) und den Schwellenkontrast 
(C) eines noch erkennbaren Details, ist als Detailinformationsgiite zu inter- 
pretieren. 

Bildwandlergerate haben aber im allgemeinen nicht ein einzelnes, sondern 
moglichst viele im Gesichtsfeld (H) befindliche Einzelheiten (h?) zu ver- 
mitteln. Unter Beriicksichtigung dieses Umstandes nimmt (3.15) mit 


ghd ( coun NS G (tein te #6)? 
pele ae rari) ce ea 


folgende Form an: 


pene yl %¢\2 D2, 
<x ag (et) i 
nt 8B : 249 
x “ni, K® ke Bia 
S:spes outs 
1+ 5 


I, ein Mass dafiir, wieviele mit einem Schwellenkontrast (C) belegte Einzel- 
heiten (A?) aus dem Gesichtsfeld (G) erkennbar sind, ist als Gesichtsfeldinfor- 
mationsgiite aufzufassen. 

Bei der Verwendung von Infrarotsichtgeraten im Anstrahlverfahren wird 
die Strahlungsdichte, unter der ein Gegenstand erscheint, durch einen Infrarot- 
scheinwerfer hervorgerufen. Unter Vernachlassigung der atmospharischen Ver- 
luste ist die Strahlungsdichte eines Gegenstandes gegeben durch 


de OF ak ne (3.18) 
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R Abstand Scheinwerfer-Gegenstand; D Durchmesser der scheinbaren Strah- 
lungsquelle (zum Beispiel Reflektordurchmesset) ; o Reflexionsvermégen des 
Gegenstandes (Lambert-Strahler) ; S;; Strahlungsdichte der Strahlungsquelle 


oe =| Ser(A) H(A) aA; (3.19) 


Ssr(4) spektrales Emissionsvermégen der Strahlungsquelle; /(A) spektrale 
Durchlassigkeit des Infrarotfilters. 

Unter der Voraussetzung ausreichender Strahlungsdichte der Scheinwerfer- 
quelle existiert bei variablem R und konstantem g ein Sichtraumgebiet, dessen 
Schirmbildleuchtdichte die Untergrundaufhellung zu vernachlassigen gestat- 
tet; mit (3.16) bzw. (3.17) gilt demnach: 


= eye at n 1 
Be are SSB (Sopa ae + Se) - (3.20) 


Aus dem Umstand, dass die Quantenausbeute der Netzhaut allein auf Kosten 
der Sehscharfe zunimmt und das Produkt OB im Leuchtdichtebereich 
10-3 [asb] < B < 108 [asb] konstant ist, geht an Bb mit (3.20) zudem die 
Forderung: 


B = 10-3 [asb] . (3.21) 


In diesem Sichtraumbereich fallt die statistische Schwankung gegentiber der 
Photonenzahl ausser Betracht. Bei Vernachlassigung der Variation des Augen- 
pupillendurchmessers mit der Leuchtdichte bestimmt sich die Informations- 
giite allein aus optischen Geratedaten: 


a? C22? = const, (3522) 
fies ae — [const , (3.23) 
G log 2 
ee a= [ (tg =~ ) ‘const . (3. 24) 


Im folgenden Sichtraumgebiet 


ne BES 25 
dar a 4 Ki. S2B(S. m, 4 K%, + Sq) (3-25) 
erfahren die Informationsgiiten infolge der mit wachsendem R abnehmenden 
Schirmleuchtdichte und der damit vermehrt ins Gewicht fallenden Unter- 
grundaufhellung Einbussen, zu deren Veranschaulichung der folgende Reich- 


weitebegriff dienen mége. 
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Verstehen wir unter R diejenige Gegenstandsweite, in der eine zu definie- 
rende Einzelheit (C2 «2) noch gerade erkennbar ist, so kann sie mit (3.15), (3.18) 
und unter Vernachlassigung der atmospharischen Verluste durch folgende 
Beziehung festgelegt werden: 


iD) @ Scr 


a Roe eee 
360602, +i m*, K sg t) Os yaaa 
( a EO ont B Ss 


les 


4. Grenzen und MOglichkeiten 


Ausreichende Strahlerleistung, Reduktion der Untergrundaufhellung sowie 
geeignete Dimensionierung der im technologischen Spielraum wahlbaren Ge- 
ratedaten sind Méglichkeiten dafiir, die in der Statistik des «Sehvorgangs» 
begriindete Begrenzung der Erkennbarkeit auf eine fundamentale Schwelle zu 
senken. Sowohl der Strahlungsdichte der Scheinwerferquelle als auch dem 
Durchmesser der scheinbaren Strahlungsquelle sind aus technologischen oder 
dem iibergeordneten Griinden Grenzen gesetzt. Den andern Méglichkeiten 
kommt daher grosse Bedeutung zu. 


4.1 Untergrundaufhellung 


Die im wesentlichen durch Streustrahlung aus dem Sichtraum sowie durch 
thermische Elektronenemission der Photokathode verursachte Untergrund- 
aufhellung kann durch geeignete Aufstellung von Sichtgerat und Scheinwerfer 
(Parallaxe) bzw. durch Kiihlung der Kathode weitgehend unterbunden werden. 

Auf diese Weise gelingt es insbesondere, den thermischen Anteil auf eine 
Schwelle zu senken, wo von dieser Seite keine Begrenzung mehr vorliegt. Die 
dazu notwendige Kathodentemperatur hangt unter anderem von der Grenz- 
wellenlange, der Schirmausbeute und der Beschleunigungsspannung ab. Bei 
Vernachlassigung der Streustrahlung aus dem Sichtraum bestimmt sich der 
durch Beseitigung der thermischen Aufhellung erzielbare Reichweitegewinn 
mit (3.26) zu: 


i) 
R, me, 4 KE > + %G 
hy = 7) 1 S : (4.1) 
m2, 4 Ke. 


R, Reichweite des Gerats mit gekiihlter Kathode; R, Reichweite des Gerats 
mit ungekiihlter Kathode. 
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4.2. Dimensionierung der frei wihlbaren Geriitedaten 


Aus (3.20), (3.21) und den sich anschliessenden Folgerungen erwachst nebst 
der Reduktion der Untergrundaufhellung die Forderung nach Strahlungs- 
dichteverstarkung; gemass (3.5) liegt sie dann vor, wenn 


Ss” Ue Lae a 1 
hk eee eg ee (4.21) 


Dadurch, dass der Kathodenempfindlichkeit (e,,), der Schirmausbeute (a,) und 
der Beschleunigungsspannung (AU) obere Grenzen gesetzt sind, kann die Ver- 
stérkung allein mit grosser effektiver Offnung des Objektivs (1/K ,) und klei- 
nem Abbildungsmafstab (m,) der Bildwandlerréhre erzielt werden. Da kleine 
AbbildungsmaBstabe eine grosse Elektronendichte und damit eine hohe ther- 
mische Untergrundaufhellung bedingen, sind auch diesem Datum, wenigstens 
bei ungekiihlter Kathode, Grenzen gesetzt (3.20). 

Im Gegensatz zu den Tagesfernrohren ist es bei Infrarotsichtgeraten nicht 
moglich, den zu beobachtenden Gegenstand mit unbewaffnetem Auge aufzu- 
suchen und ihn nach vollzogener Fixierung mit dem Sichtgerat naher zu be- 
trachten. Um dennoch eine Orientierung zu ermdéglichen, geht es mit (3.17) 
darum, aus einem méglichst grossen Gesichtsfeld (G) eine Héchstzahl von 
Detailinformationen zu erhalten. Die dadurch erwachsenden Forderungen nach 
grossem Gesichtsfeldwinkel, grosser Gesamtvergrésserung, grosser effektiver 
Offnung des Objektivs und hoher Strahlungsdichteverstirkung lassen sich 
gemass 

Dr 1 


(og 
Diges = las Ro 


A 
. 3 ; (4.22) 


25 


(U, Unscharfe des Auges in deutlicher Sehweite) allein durch hohe Bildwand- 
lerauflésung (1/U,), grosse wirksame Photokathodenflache und eine im 
Rahmen (4.21) zulassige Objektivbrennweite realisieren. Da die Fehler der 
Elektronenoptik Durchmesser und Kathodenkriimmung weitgehend bestimmen 
und die Auflésung in der Kérnigkeit des Leuchtschirms eine weitere Beschran- 
kung erfahrt, sind auch diesen Daten Grenzen gesetzt. 
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Résumé 


Le présent article a pour objet la détermination des limites et des possibilités de 
observation infra-rouge au moyen d’un dispositif optique déterminé. On décrit des 
possibilités susceptibles de réduire au maximum la limitation de l'information basée 
sur la statistique du phénoméne visuel. 


(Eingegangen: 13. Dezember 1957.) 


On Prager’s Hardening Rule 


By Ricu. THORPE SHIELD, Providence, R. I., USA), and Hans ZIEGLER, Ziirich*) 


1. Introduction 


In order to describe the behaviour of a rigid-work-hardening material, one 
needs 
(a) an initial yield condition, specifying the states of stress for which plastic 

flow first sets in; ; 

(b) a flow rule, connecting the plastic strain increment with the stress and the 
stress increment; 

(c) a hardening rule, specifying the modification of the yield condition in the 
course of plastic flow. 

It is customary to represent the yield condition as a surface in stress space, 
convex [1]*) and initially containing the origin. The current yield conditions 
for a metal are those of v. MIsEs [2] and of TreEsca [3]. The flow rule generally 
accepted is also due to v. MisEs [4]. It is justified to a certain extent by physical 
reasons [5, 1], and it states that the strain increment vector lies in the exterior 
normal of the yield surface at the stress point. As to the hardening rule, there 
are various versions in use. The rule of isotropic work-hardening [6, 7] assumes 
that the yield surface expands during plastic flow, retaining its shape and 
situation with respect to the origin. Another rule, developed by PRAGER [8], 
assumes that the yield surface is rigid but undergoes a translation in the direc- 
tion of the strain increment. This rule accounts for the Bauschinger effect 


’) Brown University, Division of Applied Mathematics. 
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_ observed in the materials in question. The main advantage of the rule is that 
_ for piecewise linear yield conditions, such as that of TRESCcA, the law exhibits a 
limited path independence of the final plastic strain with a resulting simplifi- 
cation in the mathematical analysis. 

The following sections contain a discussion of PRAGER’s hardening rule and 
its implications for special states of stress prevalent in practical applications. 

Mention should be made of the work of Hopce (see [10], for example), 
which uses a strain-hardening rule which is a combination of the PRAGER 
rule and isotropic hardening. 


2. Treatment in 9-Space 


Let us consider an element of a rigid-work-hardening solid, referred to an 
orthogonal coordinate system x;. The state of stress of this element can be 
represented by a stress point P in a 9-space o;,. In this space, the initial yield 
surface is represented by an equation 


F(o;;) = k* = const. (2.1) 


In the following, for simplicity attention will be confined to initially isotropic 
materials for which the form of the function F is invariant with respect to a 
rotation of the stress state. An initially anisotropic material can be treated in 
an analogous manner. 

The hardening rule suggested by PRAGER assumes that during plastic de- 
formation the yield surface moves in translation. After a certain amount of 
plastic flow, it is given by 

Plojn — Xx) =F, (2.2) 


where the tensor «;, represents the total translation. Because «,;, is not neces- 
sarily the isotopic tensor 6;,, where 6;; is the Kronecker delta, the material 
becomes anisotropic as a result of the hardening process. Accordingly, direction is 
important and we shall fix the coordinate system x, with respect to the element, 
small deformations being assumed. 

Due to the flow rule of v. MiskEs, the plastic strain increment de;,, con- 
sidered as a vector in the space o;,, lies in the exterior normal of the surface 
(2.2) at P. Thus, it is represented by 


“to ET ar pe (2.3) 
ik 005% 


The definition of a Prager-hardening material is completed by assuming that 
the surface (2.2) moves in the direction of de,;,; more explicitly 


da; = Cc dé; 5, ) (2.4) 
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where c is a constant characterizing the material. This work-hardening law is a 
generalization to complex states of stress of a linear work-hardening law in 
simple tension (Figure 1), which exhibits a Bauschinger effect. [The work- 
hardening modulus c, in simple tension (Figure 1) is related to the work- 
hardening modulus ¢ by c, = (3/2) c.] 


Figure 1 


Response of the material considered in simple tension or compression. 


The hardening rule described is physically acceptable because the compo- 


nents 
bin = € ErK (2:5) 


form a tensor of the second order, and the law is therefore independent of the 
particular coordinate system x; chosen. 

The scalar dA in (2.3) is determined by the condition that P remains on the 
yield surface in plastic flow. From this condition, 


OF 
(do; — do; x) ae 0, (2.6) 
and from (2.4) and (2.3) follows at once 
il (OF /00; .) do;; 
dh —— 5 See ean J uj : 
c (dF ]0c,,) (OF /0c,,) ’ (2.7) 


if the summation convention is adopted in 9-space. 
In an initially isotropic solid the yield function takes the form 


F(o;,) = G (L1(O;x), Lo(ix), I3(0;,)] , (2.8) 
where 
al 1 
I,=0;;, fee 0745 O77» Ly = > Gin Ojn On: (2.9) 


are the invariants of the stress tensor. Moreover, if the initial yield is independ- 
ent of the mean normal stress, 


Foi, + B Oix) = F(G;x) (2.10) 
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where # is an arbitrary scalar. When plastic flow has set in, the yield function 
_ becomes, on account of (2.2) and (2.8) 


Fo; — &,%) =G [L1(Oi~, — %;%), Iy(0;~ — Xe x)» I3(6i% — %ex)] - (2.11) 


From (2.10) it follows that the values of (2.11) remain unchanged when 0; 1S 
replaced by o;, +8 0;,: PRAGER’s hardening rule implies that during the whole 
hardening process yield is independent of the mean normal stress. 

- With (2.11), the flow rule (2.3) reads 


= eG sh.) 0G rer 
dex, = 4, dh ee rea +) dh. (2.12) 
Since 
or ol ol 
ae = 0ix» ie = Sigs Bigs sae (Op 4 — 049) (Oem. = pu) %, (2-13) 


we obtain from (2.12) 


0G 0G 0G 
dei, = an Oix + bir (Cin — %x) + Oly (045 — 43) (O74 — a; x)| dA. (2.14) 


Let us assume now that the physical coordinate axes originally coincide 
with the principal axes of stress. Then we have first 
Gi 0 Cah) and +o;,=0; (2.15) 
From (2.14) follows 
de,,=0 (+k). (2.16) 


I.e., since the material is isotropic at the beginning, the strain increment tensor 
is coaxial with the stress tensor. By (2.4) and (2.16), also 


da; i 0 (7 + k) . (2.17) 


The last result remains valid if the second assumption (2.15) is replaced by 


the weaker assumption 
tip =O (+R). (2.18) 


It follows that, if the principal axes of stress remain fixed in the element from the 
start, the strain increment tensor and thus the strain tensor vemain coaxial with the 
stress tensor. 

If the principal axes of stress rotate, (2.16) holds only in a first step, pro- 
vided the principal system of stress is used as the physical coordinate system. 
If (2.16) shall hold in a second step, the coordinate system must be rotated 
between the first step and the second one. This rotation, however, violates (2.18) : 
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Due to the anisotropy caused by strain hardening, the strain increment tensor is in 
general not coaxial with the stress tensor. 

Many problems of practical importance can be treated in a space of less than 
9 dimensions. In certain cases, e.g., a 3-space defined by the principal stresses is 
useful. From our last result follows, however, that this 3-space is inadequate 
where the principal axes of stress are not fixed in the element. In addition, we 
shall see in the next sections that the reduction in dimensions is not without 
influence on the form of the hardening rule. 


3. Treatment in 6-Space 


On account of the symmetry of the stress and strain tensors, the problem 
may as well be treated in 6-space. It is convenient here and particularly for the 
subsequent specializations to denote the physical coordinates by x, y, z, the 
stressés bY Op, :-., Ty2>++-, and the strains by <,,... , gp... , WHeTe temaes 
indicate cyclic permutations. 

In the new notations the yield condition (2.2) reads 


EG page Ue) Tye Bay Cpe Oe es ee (3.1) 


where T, ,, Ty, .-. have to be considered as independent variables. The flow rule 
(2.3) becomes 


Pe oe inipet te = Elie Os. a dio Wee) 
and the hardening rule (2.4) takes the form 
Akg SCASs, va05  Ckyg i Cb cmak Ring = 6 O62 ee (3.3) 
Treatment in 6-space, however, requires the elimination of the stress compo- 
nents t,,,-.., Of the strain components ¢,,,..., and of the displacements 
Cente a ia 


Because of the symmetry of the stress tensor 


ENG ys tasaehh ges tas Gg yi re) ed Oeste Seem (3.4) 


Thus, the yield surface in 6-space is given by 


Kg Oe, very Tyz — Ayz, 5) = Yala as Pesca Tyz — kyzys++,Tzy — Azy, +] =k. 
390) 
From (3.2) and (3.5) we obtain 
of 0 
deg = Gq hyn, Uys = 2idey, snd, (3.6) 
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This is the well-known result that the flow rule of v. Mises remains valid in 6- 
space, if the state of strain is represented by the engineering components &,, 
Ope 

If PRAGER’s hardening rule holds in 9-space, the yield surface (3.5) an 
_ 6-space also moves in a translation. On account of (3.3), this translation is given 
by 


eee y 


Binds, On diye = Seay, (3.7) 


im general it 1s not in the direction of the exterior normal at the point P. (3.7) is 
the form that PRAGER’s hardening rule takes in the new strain components 
in 6-space. 

It might seem that, dropping the factors 1/2 in (3.7), one might postulate 
_ the validity of PRAGER’s rule in its original form in 6-space, thereby renouncing 
its validity in this form in 9-space. Since both sides of (3.7) represent tensors, 
such a procedure would involve the sacrifice of the invariance of the rule with 
respect to rotations of the physical coordinate system. It is clear that this is 
inacceptable, and that we have to accept, conversely, the fact that the form of 
PRAGER’s rule is apt to deteriorate in a subspace. The next sections will show 
different stages of this process. 


4. Special Cases 


In many practically important cases some of the stress components are 
absent. Starting once more in 9-space, we may denote the stress components 
present by o/,, the zero ones by o;,,. The initial yield condition is then 

F(e,,, 0, = 0) = AG) =k. (4.1) 
If we are not interested in the strains ¢/, corresponding to the zero stresses o/,, 
we may treat the problem in a subspace o;,. Here, H(o;,) defines a new yield 
surface. — 

After plastic flow has set in, the yield surface is given by 

F(oj, = Hy, — Oy) = FR. (4.2) 
We will not be able, in general, to express (4.2) by means of the function H: 
in general, hardening implies a deformation of the yield surface in a subspace. 

From (4.1) follows 


; OF oH 
dé; p — 00%, dh iS 0015, i) . (4.3) 


Thus, the flow rule remains valid in any subspace. However, it supplies only the 
strain components é;, defined in this subspace although the €;,, too, may be 
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different from zero. It is clear, therefore, that even in cases where the yield surface 


undergoes a translation, it may move in a direction different from the outward nor- — 


mal at the stress point P. 


The cases where the new yield surface does not deform are those in which it 
is possible to convert the left hand side of (4.2) such that the terms —a/, 


vanish. 
When the initial yield is independent of the mean normal stress, the initial 


yield function can be written 


VOon via teeta) = Olan a) (4.4) 
where 
if oe 1 2 
fa Og Oy Og rt Ae gee | 
2 il 4 
Dee or Ge te ve 9. 0, gt Sq) ny Oe (Gey) 
if 
ia (2 Oy — Oy G;) ya ian ge Tyz Tze Try | 


are the invariants of the stress deviator. When plastic flow has set in, o,, ... , 


Tya>* have to be replaced by op. => ay, i.) Tye Gye --- =. Lt follows thar 
yield remains independent of the mean normal stress. 
By (3.6) 
Dp ee, sc8 Jy og OJ; 
den= (Seo + GE dan) Mh | ne 
Lor moe Biege Cay 
WY ys is (a, cae yi OTs oe, ee poe 


where, on account of (4.5) and under the assumption that plastic flow has taken 
place, 


0 Jo il 
jak = F(Z (Oe — a4) — (04 — 04) — (0p — 24) one, 
=2 (ge Kya ) 
OSs t ( 2 2 9 
00, ae) [2 {On — Or) a (o, > Oy) — (o, a G,)? ar 4 (o, sq Oy) (o, ,) 
— 2 (0, — a) (Gy — tg) — 2 (0, — a4) (0y — &)] oe 
1 
ae S [2 (le eal ar hea} (Vo Ricg! he ’ 
OJ; Pi 
Ft =~ = (a ~ %) — (04 a4) — ee 


=e 
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a 
rR 
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From (2.10) follows 
Sp Oa OPE G Tyg 20g grnen) =H (Og. Oy,» oye gi Oyegsses)» (4.8) 


_ Hence, in cases where one of the normal stresses, e.g. O,, 1s absent, the term 
— «, can be eliminated by addition of «, to the normal stresses. 
If t,, = T~ = 0, (4.6) and (4.7) yield, in connection with (3.7), 


QVyz=4y,2=0 (Yy,=Ver= 0). (4.9) 


This is the proof, by complete induction, that 


1 1 
bye = 7 oVyz = 9, tre = 7 oVee = 9 Gigebte Re Be E10) lee (4.10) 


i.e. that, in the absence of at least two shear stresses, the corresponding shear 
strains and thus the corresponding displacements of the yield surface, are 
_ absent. This result might have been inferred from the symmetry (Figure 2) of 


Figure 2 
State of stress with T,, =Tz, = 0. 


the state of stress with respect to the middle plane x, y of an element. It is clear 
that a similar result does not hold in general if only one shear stress is absent; 
on the other hand, (4.10) is the reason why PRAGER’s hardening rule applies 
without change of form in the space of principal stresses. 

While the foregoing results hold for any form of the yield function (4.4), 
we can obtain some more results by restricting ourselves to the more common 
types of g. In v. MiskEs’ case, the function g reduces to J,, and the initial yield 
condition is 

h=7%, (4.11) 


where 0, is the initial yield limit in simple tension or compression. In TREsca’s 
case, g also depends on J3. Here, the yield condition is better discussed in terms 


of maximum shear stress. 
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In (4.10), it is essential that two shear stresses vanish: in cases where one — 
shear stress only, i.e. t,,, is absent, it may be impossible to eliminate the cor- — 
responding displacement, «,,, from the right hand side of (4.8), and this means — 
that the yield surface deforms in the process of hardening. In TRESCA’s case, — 
this actually happens. However, if v. MisEs’ yield condition applies, it follows — 
from (4.6), (4.7) and (3.7) that 


1 
Cny = = oVay = 9 (Try =0). (4.12) 


If o, = o, = 0, we obtain from the same relations 
Eye e,, Ady — Oe 1G, 0, = Op. (4.13) 


provided the material obeys the yield condition of v. Misss. It is easy to see, 
however, that under TREsca’s yield condition (4.13) does not apply. In the last 
case it is not possible to eliminate both «, and «, from the right hand side of 
(4.8): the yield surface deforms in the course of hardening. + 

In the next sections we shall make use of the fact that the plastic volume 
change of the material is zero. That is, 


été t+ée,=0, (4.14) 


a consequence of the flow rule and the independence of the yield function on the 
mean normal stress. 


5. Plane Strain 


Here, per definitionem, 


Toa = tage 0, 2, = 0. (5.1) 
From (4.10) we obtain 


Vue = Vig = 0. (5.2) 
It follows from (3.7) that the yield function has the form 
5 
(6, — 6 &g, Fy — CEy, Oz, Tay — ZC Yes) i (5.3) 


On account of (3.6), 


and by (4.14) 
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Thus, we finally get the yield condition 
hon — ces, O, + C Ez, ite C yey) == RA. (5.6) 


_ It follows that the yield surface moves in a translation, but, on account of the 
factor 1/2 in the last argument, not in the direction of the exterior normal. 

In many cases where the original form of PRAGER’s hardening rule does not 
hold in the subspace appropriate to the problem, this defect can be remedied 
by a transformation. Here, for instance, the new quantities 


s 1 
ig = V2 Try » SPIES Y (5.7) 


_ can be introduced. The yield condition (5.6) becomes 
D(z — C Eg, Sy + CEg, tay — C Bey) = R*, (5.8) 


and the form of PRAGER’s rule in 9-space is regained. 
In order to specialize (5.8) for v. MisEs’ yield condition, we start from 
(4.5) and obtain first 


1 

f=h= 5 (o, — ez)" + ——- S. (= C Ey) (o, — ¢ &,) Sai di 

(5.9) 
1 2 
oe (tas — = Yay) yg 
Carrying out the steps (5.3) through (5.8) with (5.9), we get 
“i 

= [(Ox Sale Ex) e (Gy si C Ey) |? +2 (toy gue Bey)” =) ee 0 - (5.10) 


Figure 3 shows the yield surface, a circular cylinder of radius (2/3)"” a4 the axis 
of which is parallel to the plane o,, o, and bisects the angle between the axes 


0, and oy. Os 
In order to specialize (5.8) for TRESCA’s case, we start from the principal 


stresses 


Cn Oe 2 1/2 
O42 = ae + |( a x) + 72,| Og = 0, (Sef) 


where a; lies between o, and o,. The material yields initially when the maximum 
shear stress reaches a critical value, 


Bi os eas an OF 5.12 
(ahold Sp (3.12) 


After plastic flow has set in, (5.12) takes the form (5.10) with of replacing 


/ 
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4 o2/3 on the right hand side. The yield surface is a circular cylinder (Figure 3) — 
of radius oo|V2, moving in the direction of the outward normal. 

Subcase a: If t,, = 0, (4.10) yields y,,, = 0. The problem can be treated in a — 
(oz, 5,)-plane and is a degenerate case of principal stress space with the prin- 
cipal strain ¢, zero. The yield locus is the strip obtained by bisecting the cylin- 
der of Figure 3 parallel to the plane o,, o,. 


Figure 3 


Yield surface in plane strain. 


Subcase b: If o, = 0, the problem can be treated in a (o,, T,,)-plane. The 
yield locus, obtained by intersecting the cylinder of Figure 3 with a plane paral- 
lel to the plane o,, ¢,,,, does not move in the direction of the exterior normal. 
However, if besides (5.7) the transformation 


1 a 
a= 5 Om ep =) 2 &e (5.13) 
is used, the yield condition (5.10) with o, = 0 becomes 
Bs 
(Se eee Cy)” ae (toy ae ay)” Ries 9 - (5.14) 


The yield locus is a circle of radius (2/3)'/* og in v. MisEs’ case and (1/2)"? a, in 
TRESCA’s case, moving in the direction of the exterior normal. 


6. Plane Stress 
Here, per definitionem, 
Op = Vpn ST gp SO (6.1) 
On account of (4.10), 
VYuz =Ver = 9, (6. 2) 
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_ and from (4.14) follows 
& = —(e,+ &,). (6.3) 
Hence, the yield function is 


ak 
g|oe— cee Oy — C Ey, C (Ey + &), Tea pou (6.4) 


_ After subtraction of c(e, + ¢,) from the normal stresses in accordance with 
(4.8) we obtain the yield condition 


1 
ilo. Sant s,), 0,— 6 (&, + 283), Tey — 5 c Yey| = Fe tO 5) 
Again, the yield surface moves in a translation, but not in the direction of the 


- exterior normal. 
If we make use of the transformations 


1 = 

Ss => (G2 + Oy), See le Gels te, =\3 Tey 
(6.6) 

1 
Cg = Ey + Ey, ula (Ey — &,) , Sen 3 Vays 
the yield surface becomes 
3 3 3 

b(s,-— = ces, Sree ae ten — © Ben) = (6.7) 


and moves in the direction of the outward normal. 
In the case of v. Mises’ yield condition, (6.7) takes the form 


p= (s, — = A 28 (s, — = wal + (t,, — a EAN =Gi« (0.9) 


The yield surface (an ellipsoid in the original stresses) is a sphere of radius op. 
In TREsSCA’s case, we obtain three sheets 


gs? + 72 \1/2 
Ae 2 | =O 


s? + 72 1/2 
ee ee es 


3 


Thus, the yield surface (Figure 4) is a circular cylinder, closed by two circular 
cones such that the intersections of the whole surface with the middle planes 
parallel to se, s, and sg, ¢s, are regular hexagons with sides oy. 

It often happens that the yield surface contains singularities. In 3-space 
such singularities are edges where two smooth surfaces intersect (in Figure 4 
the circles of intersection between cylinder and cones) or corners where more 
than two smooth surfaces meet. An isolated singularity or vertex (in Figure 4) 
the vertex of either cone) may be considered as a limiting case of a corner. 
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If P is to remain in a corner or vertex, the yield surface must move with P; 
hence, the displacement d« coincides with the stress increment do. The strain 
increment de follows from da, since the components of a are the factors of c in 
the general equation [here (6.5) or (6.7)] of the yield surface. If dx = c dé, then 
P remains in the corner or vertex as long as do lies in the pyramid or cone en- 
closed by the exterior normals of the yield surface in the vicinity of the singu- 


Figure 4 Figure 5 
Tresca yield surface in plane stress. v. Mises yield locus in plane stress with Ty, = 0. 


larity. If dx + c de, the pyramid or cone is obtained from the vectors dx instead 
of the normals. 

If P is to remain in an edge, the components of do and dx normal to the 
edge must coincide. Further, according to the flow rule, dé has no component 
along the tangent of the edge. Thus, two components of dx and the third one 
of dé are known as soon as do is given. The remaining components follow again 
from the equation of the yield surface. If dx = c de, P remains in the edge as 
long as do lies in the wedge enclosed by the normals of the yield surface in the 
vicinity of the edge. If dx + c de, the vectors da take the function of the normals. 

Subcase a: If t,, = 0, (4.10) yields y,,,= 0. The problem can be treated in 
a (0,, 6,)-plane, the section of principal stress space by the plane o, = 0. 

In v. MisEs’ case, the yield locus is the well-known ellipse illustrated in 
Figure 5 with the equation 


[Oz cai (2 Ex a &y) Sie [Oy Pei (& a 2 ey? | 
; (6.10) 
—[o,—¢ Zez te) [0,—¢les+2e)—=o8. | 


Comparing the vectors dx,, dx, and de,, de, following from (6.10), we easily 
obtain 


do, = 3¢(t,—A,) dA, doa, =3c¢(o,—4a,) dd. (6.11) 


Hence, the ellipse moves in the direction of the radius C P. 
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Using the transformations (6.6), we obtain a representation in a plane 

Sz, S,. Here, the yield locus is a circle of radius o,, moving in the direction of the 
outward normal. , 

. In TREsca’s case, the yield locus in the (¢,, 0,)-plane is the hexagon of 

Figure 6. The arrows indicate the vectors dx, of constant direction on each side, 

_ given by the line connecting the center of the side in question with the center 


Figure 6 Figure 7 
Tresca yield locus in plane stress with T,, = 0. Yield locus of Figure 6 in sg, Sy. 


_of the opposite side. If do lies in one of the shaded regions, the hexagon moves 
with P. Taking dx,, dx, from (6.5), solving for de,, de, and setting dx = do, we 
obtain 


de, =~ (2do,—do,), dey = =~ (2 do, — doy) (6.12) 
_ for the strain increment in a corner of the hexagon. 

In the (sg, s,)-plane the yield locus becomes a regular hexagon (Figure 7) 
with side o, moving in the direction of the exterior normal. It is clear that this 
hexagon is the section of the yield surface of Figure 4 with the plane of sym- 
metry parallel to sz, s,,. 

If P remains in a corner, 

d= 7 dy, de, = =~ ds,; (6.13) 
thus, de has the same direction as ds = da. 

Subcase a: If o, =0, the problem can be treated in a (6, T,,)-plane. How- 
ever, since (4.13) only applies in v. Misks’ case, the yield locus deforms in any 
other case in the process of hardening. 

In Tresca’s case, for instance, the initial yield locus is the intersection of 
the yield surface of Figure 4 with one of the vertical planes passing through C 
and touching the two circular edges. The displacement of C, however, generally 
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does not lie in this plane, and it becomes clear, therefore, that the yield locus 
deforms. 

In v. Mises’ case, (4.13) applies. Hence, the yield locus does not deform. 
This follows also from the fact that the yield surface corresponding to Figure 4 
is a sphere. The yield locus is the circle 


(o. = c ez) ae (t., = c Ben) = On (6.14) 


in the plane o,, ¢:,, where te, and gz, follow from t,, and y,, by means of (6.6). 
The circle (6.14) is of radius og and has been discussed by PRAGER [9]. It moves 
in the direction of the outward normal. 


7. Another Special Case 


In certain cases, e.g., if a cylinder is subjected to torsion and simple tension, 
we have 0, = dy = Try = 0. From (4.12) through (4.14) follows 


a 
by = by = a Ee Joy =U; (7.1) 


provided the material obeys v. MiskEs’ yield condition. For any other yield 
condition (4.12) and (4.13) do not generally hold; hence, the yield surface in 
3-space deforms in the process of hardening. 

Restricting ourselves to materials obeying v. MIsEs’ condition, we start 
from (4.5), i.e., from the yield condition 


fae Gp mic ee 


Inserting (7.1) and subtracting c ¢,/2 from the normal stresses, we obtain in this 
particular case 


Sa (c. cs = ce.) ote t3 = = Cus) +3 (Res Ta 2 Yea). = o% . (7.3) 


This yield surface (an ellipsoid of rotation) moves in a translation, but not in 
the direction of the outward normal. 


If we use the transformations 


PN. ay iy 1 
bys = [3 Tye s tee =3 Tz a» ott hin fav ted PRT Te AEA (7.4) 
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the yield surface becomes a sphere 


2 


p= (0,— 3 ce.) + (tz eee) + (tee 5 Bes) = 08 (7.5) 


of radius o9, moving in the direction of the exterior normal. 
Subcase a: Ifo, = 0, (4.13) and (4.14) yield ¢, = 0. For v. Msgs’ yield condi- 
tion, the problem can be treated in a (t, ,, T,,)-plane. The yield locus is the circle 


2 il, 2 


(tye ere) (t2— Fo Yes] = Op (7.6) 


of radius 3~'/? 6, moving in the direction of its outward normal. 
, It is interesting to note that under TrEsca’s yield condition the yield sur- 
face deforms even in this comparatively simple case and that y,, does not 
remain zero. 
Subcase b: If ty, =0, (4.10) yields y,, = 0. If v. Msgs’ yield condition 
holds, the problem can be treated in a (a, ¢,,,)-plane. The yield locus is the circle 


2 


(c.-+ce,) + (te2— 3 682) = 95 (7.7) 


of radius o), moving in the direction of the exterior normal. 
It is clear that the last result, apart from the difference in notation, is the 
one already obtained in the last subcase of section 6. 


8. Conclusion 


In sections 5 through 7 we have encountered the various possibilities discus- 
sed already in section 4. In most of the cases considered the yield surface moves 
in a translation, and a simple transformation at most suffices to make the 
original form of PRAGER’s rule apply in the subspace appropriate to the problem. 
In certain cases, however, the yield surface deforms in the course of hardening. 
Incidentally, these exceptions occur, as far as our examples are concerned, in 
those cases where at the same time (a) more than one normal stress is different 
from zeroand (b) the material obeys TRESCa’s yield condition. 

In section 1 the advantages of TrEsca’s yield condition have been emphasiz- 
ed. The deformation of the yield surface which has been found in many cases 
represents a serious drawback in this respect. 
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Zusammenfassung 


Um das Verhalten eines Metalls mit dem Spannungs-Dehnungs-Diagramm 
der Figur 1 unter einem beliebigen raumlichen Spannungszustand zu beschreiben, 
wird neben der Annahme isotroper Verfestigung, welche aber den Bauschinger- 
Effekt nicht erklart, die Pragersche Verfestigungsregel [8] verwendet. Es wird 
hier untersucht, welche Formen diese Regel in den wichtigsten Spannungsraumen 
von weniger als neun Dimensionen annimmt. 
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On the Free Convection from a Horizontal Plate 


By KEITH STEWaRTsON, Durham, England?) 


The free convection of heat from a heated vertical plate in a fluid has been 
extensively studied in recent years. A review of the work done has been given by 
Squire [4]?) and subsequently numerical solutions of the governing equations 
has been given by Ostracu [2] for a wide range of values of the Prandtl number o. 
The convection takes place in boundary layers originating at the lower edge of 
the plate. Fluid is drawn into them, is heated and gaining buoyancy moves 
upwards. On the other hand if the plate is cooled relative to the surrounding fluid 
the situation is reversed for the boundary layers originate at the top of the plate, 
and the fluid drawn into them is forced downwards. When the plate is inclined 
to the vertical there is no change in the flow pattern, since the vertical buoyancy 
force has a component along the plate which drives the fluid thus generating the 
boundary layer. However, if the plate is horizontal the buoyancy has no compo- 
nent along its length and the boundary layer, if it exists, must be of a different 
character. 

’) Department of Mathematics, The University. 

“) Numbers in brackets refer to References, page 281. 
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In this paper we shall examine the boundary layer on a horizontal plate 
Showing that if the plate is heated and faces downwards or is cooled and faces 
upwards a solution can be found which is in moderately good agreement with 
experiment. On the other hand if the plate is heated and faces upwards or is 
_ cooled and faces downwards the boundary layer problem is not properly posed. 

The appropriate equations for a heated plate inclined at an angle a to the 
horizontal are [4] 


du Ou Op : 0?u 
eu Ox + Qu oy aad a OX ee ine (1) 
0 
0=— 5°. —egcosa, (2) 
Ou ov 
mee Teak t (3) 


in which x measures distance along the plate and upwards, y distance perpendi- 
cular to the plate and upwards, u,v are the components of velocity in these 
directions; p the pressure, g@ the density, T the absolute temperature and wu the 
coefficient of viscosity. In deriving these equations it is assumed that both the 
viscosity and the difference between the ambient temperature T, and the plate 
temperature T, are sufficiently small. Let 


P = Po — On V COSK — Oog XSINA + Oy P, (5) 


where gy, is the constant density of the fluid outside the boundary layer and 
is a constant. Since the boundary layer is thin, it follows from (2) that oP, 
which is the pressure increment inside it, is small. From the equation of state [4] 


we also have 
@—Q=—Bea(IT—-T), (6) 


where f is a constant. In particular for a gas B7,=1. The first two equations 
of motion now reduce to 


! Or 02u 
Ou Ou (7) 


sce Felina 5 eet; rie Pgisith a. (Liem Ba) iPass 


where » is the kinematic viscosity. 

The boundary conditions are that w=v=0, d= Datei On and 470) 
T=T,, P=0 outside the boundary layer or, effectively, as y > oo. By them- 
selves these are not sufficient to solve the problem completely since the equations 
are parabolic and initial profiles of T and wu must be prescribed as well. If aw +0 
the difficulty may be overcome as follows. From (7) w must be positive if the 
plate is heated and therefore the boundary layer must originate at the lower 
edge of the plate where it must contain a zero mass of fluid. Similarly if the plate 
is cooled it must originate with zero mass at the upper edge of the plate. 
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If ~ = 0 however it will appear below that the boundary layer does not always 
exist and even when it does it may originate at any value of + growing either 
in the direction of increasing x or in the direction of decreasing x. We shall 
assume the most obvious possibility, that the boundary layers originate at each 
of the two edges of the plate and grow towards the centre. ScHmiptT [3] has 
published a photograph of a heated horizontal plate in which the boundary layer 
is thicker at the centre than at the edges. It may be regarded as supporting but 
not conclusive evidence for the present assumption. 

The solution of the governing equations for a horizontal plate may now be 
attempted. Take the origin of x to be one of the edgés and write 


[Leeds 1/5 
ym LAL Te TA PERE wa (489° B g cose |T,— Tol} F(a), 


x vy (9) 
T=T,—|Ty—T | 6(y) sgny, P= —[» B?g? x cos?a (Ty — Ty)*]** G(n) seny, 
where y is the stream function defined by 
_ oY 2 ip 
ee Oy eee ae Y: 
so that 
u = [x v B® g? cos? (T,— To)?}° F’(n), 
(10) 


1 
v= bess [x8 v8 B g cosa |T,—T)|]1/5 (F—7F’). 


The factor cos«, equal to unity when the plate is horizontal, has been retained 
for convenience when inclined plates are discussed below. The equations for F, G, 
0 are, since « = 0, 


vi 3 "” 1 / 2 z / " 3 / 
Hi, FF F822 (6 9G), “G'=0, 6’ p— o Pe Sore 


with boundary conditions 


6= +sgny, Na TON Ee pie Oe 
| aa 


06+0, P+>0, F’+>0 as 7>o. 


The special form of these equations shows that there are only two distinct 
cases according as 6 = +1, or = —1 at 7n=0. 

First suppose that = —1 at 7 =0, corresponding to a heated plate facing 
upwards or a cooled plate facing downwards. In this case (11) has no solution 
for on integrating the equation for F twice 


| 
, 3 2 i 
Fn) + yy Fen) = — 3 f an [ (F246 96) an. 
Hence , ‘ 
CO lee) me.) 
SSP fa (io) We ee vy 19 j 2 
= == | dy | (2F4G—nG)dn= ai dn (2 F’2+ 420’) dn, (13) 


0 n 0 
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which is a contradiction, since from the equation for 6, 6’>0. Thus the problem 
_ of a heated horizontal plate facing upwards has not been properly posed and most 
probably the plate does not have a boundary layer. 
Second suppose that 9= +1 at 7=0, corresponding to a heated plate facing 
_ downwards or a cooled plate facing upwards. Then (12) presents no difficulty 
and has been integrated numerically taking o = 0-72, its value in air, it being 
found that 


9’(0) = —0-358,  F"(0) = 0-971, G(0) = —1-73,  F(co) = 2.30. (14) 


The form of the boundary layer is slightly different from that of a heated 
vertical plate because the buoyancy acts indirectly. It induces a transverse 
pressure gradient whose variation along the plate drives the fluid. When 6(0) = —1 
the pressure gradient along the plate which is induced is adverse thus explaining 
why no solution can be found. It is to be expected, and confirmed by (9), that 
since the buoyancy acts indirectly the boundary layer is weaker and thicker than 
it is on a vertical plate. 

The flow pattern is similar to that arising from a finite circular cylinder rotat- 
ing about its axis in a fluid which is rotating about the same axis but with a 
different angular velocity. Some problems of this kind have already been ex- 
amined [5] and it may be shown that the basic equations are the same if o=1 and 
the two angular velocities are nearly equal. Indeed the present paper was stimu- 
lated by a discussion on rotating fluids with Professor Squire. During it he 
showed me Scumipt’s photograph of a heated horizontal plate pointing out that 
just as there is a fundamental difference between the flow inside and outside a 
rotating cylinder, so is there a difference between the flow above and below the 

late. 

The solution given by (9) and (14) refers to a boundary layer originating at 
one edge of the plate and there will also be another equal boundary layer origi- 
nating at the other edge. The fluid in each of these layers is moving towards the 
centre where a collision between them must occur. Since the fluid in each layer 
has momentum, fluid must then be expelled moving away from the plate in a 
- two-dimensional jet. The interaction between the boundary layers is similar to 
that in the theory of a rotating sphere already discussed [5]. Beyond a distance 
O(v%/5) from the centre, the two boundary layers are independent. Within this 
distance each boundary layer divides into two layers, of which the outer layer is 
a jet, being controlled solely by the inertia terms and the inner of thickness 
O(yv%/5) is a conventional boundary layer with a mainstream and must separate. 
If, as is believed, the effect of the separation is small the inner layer may be 
neglected, when the initial profile of the jet is the same as the profiles of the two 
boundary layers just outside the interaction region when taken together. 

The total heat transfer from unit length of a plate of breadth / is 


1/2 
(oy ov y=0 
0 


Substituting from (9) and (14) we have 
Le ae ed a (T1—Th) 6 y= 0.841 (Ra)¥5, (15) 


Nex = : 
bee, (1.15) \ ” 


where Nu, Ra are the Nusselt and Rayleigh numbers respectively.. 
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FISHENDEN and SAUNDERS [1] have reported experiments on horizontal plates. — 
These showed considerable qualitative differences according as the plates were — 
heated or cooled, faced upwards or downwards, in agreement with the predictions ~ 
of this paper. For heated plates facing downwards, when a quantitative compari-— 
son can be made with the results of this paper, they concluded that the flow is 
laminar throughout the range of their experiments and that 


Nu = 0-25 (Ra)"4, (16) 


According to the graph of Nu against Ra which they display however the slope 
of the curve is almost exactly 1/5 in 10®°<Ra<108 thereafter increasing to 
1/4 as Ra increases to 101°. It may be that at the upper end of the range of Ra the 
flow although not turbulent was slightly unstable. Another possible explanation 
of the increase in the index is that the temperature difference 7,—7T, was 
ultimately over 500°C when the effects of the variation of viscosity with temper- 
ature and of dissipation may be important. Assuming that Nu «(Ra)/* throughout 
the closest fit with the curve is obtained when 


Nu = 0-81 (Ra). (17) 


The agreement is much better than when (16) is used. Although (15) and (17) 
agree quite well it must be remembered that the experiments were done with 
plates which were nearly square, while (15) is appropriate only to an infinitely 
long strip. The theoretical result for a square must have the same index as (15), 
from dimensional considerations, but since the square has four edges it would 
have a slightly higher coefficient. 

Let us suppose now that the plate is inclined to vertical. The argument of 
this paper may still be used if 


ORR . , 
aa >B eg (T—T)) sing, 
1. €., 
[v B? g2 cos?@ (T, — Ty) 22/5 4—3/5 G >B g|T,—T)| 6 sine. (18) 


Thus no matter how small « is then ultimately the component of the buoyancy 
along the plate must dominate the induced pressure gradient. To find when (18) 
holds let us assume that G = O(@) and that S means more than 10 times as large. 
The condition is then that 

Eel | cos*g ae 

l a | 105sin'e Raf * 


Thus the present theory is valid everywhere on the plate when Ra = 105 provided 
that the inclination of the plate to the horizontal is less than one degree. 

On the other hand the condition that the component of buoyancy perpen- 
dicular to the plate be negligible can also be estimated. It is that 


[v B® g? cos? (T,—Ty)?]*5 x-25 G <B g|T,—Ty| sing. (19) 
It is clear from (19) that for any « </2 the component of the buoyancy perpen- 
dicular to the plate must dominate the component parallel to the plate sufficiently 


near to the origin of the boundary layer. If we require (19) to hold over OE os: 
the plate then 


108 cos4g < sin’g Ra, 
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whence if Ra=10°, » >70° approximately. Thus for plates inclined to the 
vertical the transverse pressure gradient induced by the buoyancy is important 
near the origin of the boundary layer over a large range of values of the Rayleigh 
number. Accordingly the boundary layers on either side of sucha plate should be 
slightly different near their common origin. 

The solution given here is valid if three assumptions are justified. First the 
boundary layers must be centred on the edges of the plate. Second the assumption 
implicit in (1)—(4), namely that v <u must hold which means, from (10) and (15), 
that 


ia tae (Ra) a bles (Ra)5, i.e. (7) S> (Ra)-¥8, (20) 


Third the convection currents generated by other heated bodies in the neigh- 
bourhood must be negligible when compared with the velocities in the boundary 
layers. If the heated plate forms the ceiling of a closed room, of which the walls 
and floor are maintained at the same temperature, the only convection currents 
arise from the ceiling and do not modify the solution given here provided that (20) 
is satisfied. If however it is a face of a heated rectangular bar convection currents 
are set up from the vertical sides which are larger than the current set up from 
the horizontal face. There is a possibility therefore that in certain circumstances 
they may not be negligible in comparison with the velocities in the boundary 
layer. It is known [4] that the convection currents induced by a heated vertical 
plate are of order v/~—1 (Ra)'/4 assuming that a typical length in the vertical 
direction is also /. Further at a lower corner of the bar the convection currents 
outside the boundary layer may be singular the velocities being at worst of order 
v3, where 7 is the distance from the corner. Hence the velocity of the convection 
currents just: outside the boundary layer are 


y (1\I1/ 1 
O [+ (=) (Ra) 
and are small compared with u if 


a 


ike (F)" (Ra)? SS + (=)" (Ra)"4, i.e. ; S> (Ra) — 982, (21) 


PONG 


The difference between the indices in (20) and (21) is sufficiently small that the 
effect of these external convection currents may be neglected as soon as the 


boundary layer approximation is valid. 
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Zusammenfassung 


Die freie Konvektion an einer beheizten horizontalen Platte wird nach den 
iiblichen Methoden der Grenzschichttheorie untersucht. Weist die warmeiiber- 
tragende Plattenseite nach oben, so ist das Problem auf dieser Basis unlésbar, 
wohl aber findet man eine widerspruchsfreie Lésung, wenn diese Plattenseite 
nach unten weist. In diesem ist die Nusselt Zahl, die den Warmeiibergang kenn- 
zeichnet, proportional der 1/5-ten Potenz der Rayleigh Zahl, welche die Beheizung 
charakterisiert. Die Ubereinstimmung mit dem Experiment ist befriedigend. 


(Received: February 27, 1958.) 
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Friihjahrstagung der Schweizerischen Physikalischen Gesellschaft 
vom 3. bis 4. Mai 1958 in Romont 


Berichte iiber die Sitzungen fiir angewandte Physik und Mathematik 


Quelques considérations sur le développement des compresseurs 
axiaux transsoniques. Par P. Scowaar, Dammarie-les-Lys, France’). 

Pendant longtemps, les constructeurs de compresseurs axiaux se sont systé- 
matiquement abstenus de dessiner leurs aubages pour des nombres de Macu 
supérieurs 4 ceux reconnus comme limites sur la base des essais de grilles planes?). 
Ces limites, au-dessus desquelles on fait apparaitre une chute rapide du rendement 
de compression des grilles, se situent approximativement dans un domaine sub- 
sonique 0,7 < My. < 0,9 des vitesses relatives de l’écoulement d’amont, les 
valeurs inférieure et supérieure correspondant respectivement a des sections 
d’aubage voisines du moyeu et de la périphérie d’une roue mobile, ceci pour des 
rapports diamétraux de l’ordre de 0,4—0,5, c’est-a-dire pour les premiers étages 
de compresseurs couramment utilisés dans les turbomachines d’aviation. 

L’évolution de la technique de ces machines vers la réalisation de débits 
spécifiques de plus en plus élevés, s’est traduite par l’emploi généralisé d’écou- 
lements giratoires rotationnels A l’entrée des compresseurs. Ces types d’écoule- 
ment se distinguent des écoulements irrotationnels classiques par une diminution 
de la composante axiale de la vitesse et une augmentation du moment cinétique 
du rayon intérieur au rayon extérieur de la veine. Ils permettent ainsi d’éviter 
les vitesses relatives trop élevées au sommet des ailettes et de réaliser des débits 
spécifiques optima par une utilisation intégrale des limites subsoniques des 
nombres de Macu sur toute la hauteur des ailettes. Il en résulte qu'on ne peut 
obtenir une augmentation appréciable du débit spécifique qu’en dépassant ces 
limites. Un gain de l’ordre de 20% peut étre réalisé en conservant les conditions 


) Société Nationale d’Etude et de Construction de Moteurs d’ Aviation. 
2 z soe os ‘ . : 
) Nous ne considérons pas ici les compresseurs axiauX supersoniques au sens habituel du 


terme, dans lesquels on tente de mettre en jeu des nombres de Macu nettement supersoniques sur 
toute la hauteur des ailettes d’une couronne fixe ou mobile. 
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d’écoulement subsoniques classiques au moyeu et en autorisant un nombre de 


_ Mach relatif légérement supersonique My, =1,1—1,3) au rayon extérieur de 


Vaubage mobile. Ceci revient simplement a rétablir lécoulement giratoire irrota- 
tionnel a l’entrée de la roue, voire a supprimer le distributeur d’entrée, donc a 
maintenir la vitesse débitante constante sur le rayon. 

Cest a cette définition que répond notre conception initiale de l’étage trans- 
sonique a taux de compression modéré. Un relévement du niveau de Macu relatif 
au rayon intérieur permet un gain supplémentaire de débit et une augmentation 
proportionnée du taux de compression. 

Le caractére essentiel de l’étage transsonique ainsi défini consiste dans la 
présence d’ondes de choc détachées du bord d’attaque des ailettes mobiles sur 
une partie plus ou moins importante de leur hauteur et, pour les taux de com- 
pression élevés, dans l’utilisation des chocs pour produire une partie au moins du 
travail de compression de l’étage. 

Quelques propriétés importantes de l’écoulement 4 ondes de choc détachées 
peuvent étre déduites d’une analyse qualitative du fonctionnement transsonique 
de la grille plane. Dans l’hypothése de conditions d’écoulement périodiques le 
long du front de grille, on montre que le passage transsonique du nombre de Macu 
d’amont M,, s’effectue de maniére continue si la section de passage minimum de 
la grille ne bloque pas. La croissance de M,, en supersonique peut étre limitée par 
deux circonstances. Une limite inférieure M, may, COrrespond a l’apparition d’une 


configuration de choc en fourche dans la grille. Une limite supérieure WV mar? qui 


ne peut en aucun cas étre dépassée dans le schéma de fonctionnement de la 
grille (abaissement progressif de la pression statique d’aval) qui correspond 
aux conditions usuelles d’aspiration d’un étage, est donnée par la condition 
M xomax Sinn, 

L’évolution transsonique de l’écoulement est profondément influencée par la 
présence des couches limites, dont l’interaction avec les chocs prend une impor- 
tance particuliére dans la grille, non seulement par les décollements et les pertes 
qui peuvent en résulter, mais aussi par l’effet de striction supplémentaire que 
lépaississement consécutif de la couche limite d’extrados provoque au droit des 
cols de la grille. 

Il y a lieu de compléter les considérations sur le fonctionnement transsonique 
de la grille plane par quelques indications sur les influences tridimensionnelles, 
qui jouent un réle important dans la grille circulaire. Le passage transsonique 
s’effectue plus aisément dans un aubage de compresseur que dans une grille 
plane, car les conditions d’écoulement les plus critiques ne s’y produisent pas 
simultanément, comme dans cette derniére, sur toute la hauteur des ailettes. De 
plus, les écoulements giratoires irrotationnels engendrent des déplacements 
radiaux des lignes de courant vers l’extérieur de la veine a la traversée de l’aubage 
mobile, qui favorisent le déblocage des sections voisines du moyeu et tendent a 
abaisser le niveau de Macu au droit des chocs, tant a la périphérie qu’au rayon 
intérieur de la roue, réduisant ainsi l’importance des décollements et des pertes 
dues a l’interaction chocs—couches limites. 

Finalement, il importe de choisir une évolution aéro-thermodynamique glo- 
bale du fluide telle qu’il en résulte un dessin d’aubage transsonique a cambrure 
modérée. Des profils 4 faible épaisseur relative, avec maitre-couple reporté loin en 
arriére du bord d’attaque, et l’adoption de charges aérodynamiques modérées 
contribuent efficacement A éviter les grandes survitesses d’extrados. Ces condi- 
tions ne sont pas toujours compatibles avec le but final recherché, de sorte que 
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le dessin d’un étage transsonique procéde généralement d’un compromis qui 
affecte plus ou moins son rendement. 

Les données suivantes ont été adoptées pour le dessin d’un des étages trans- 
soniques expérimentaux réalisés 4 la SNECMA: 


My,= 9,700, My, = 1,324, vitesse périphérique U,=422 m/s, 
rapport diamétral y= 0,45, taux de compression 7, = 1,26, 


débit rapporté a la section frontale D,= 160 kg/m?s (INA O). 


Aux essais, le débit et le taux de compression prévus ont été réalisés au 
régime d’adaptation avec un rendement adiabatique de 84% de l’ensemble (distri- 
buteur amont + roue mobile + redresseur) et une marge de pression au pom- 
page largement suffisante (my = 1,385). L’analyse détaillée des conditions 
d’écoulement dans cet étage nous a permis de confirmer dans une large mesure 
la validité des considérations présentées dans cette bréve étude. 

Une publication détaillée est prévue dans «Technique et Science Aéronau- 
tiques». 


Uber Reibung und Grenzschicht-Stabilitat an rotierenden angeblase- 
nen Zylindern. Von J. ACKERET, Ziirich®). 

In Zusammenhang mit Untersuchungen tiber Sekundarstro6mungen haben 
wir (Verfasser und H. SPRENGER) einen einfachen Versuch gemacht, der auch 
selbstandiges Interesse bietet. 


0 05 10 15 20 25 30 


Reibungsdrehmoment am radial angeblasenen rotierenden Zylinder. 


Ein rotierender Zylinder (siehe Figur) wird radial am ganzen Umfang ange- 


blasen und das Drehmoment bei verschiedenen Anblas- und Drehgeschwindig- 
keiten gemessen. 


8) Institut fiir Aerodynamik der ETH. 


7 
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Im laminaren Fall folgt die Grenzschicht einem exakten und einfachen 


_ Gesetz: sie ist auf der ganzen Zylinderoberflache von konstanter Dicke, falls der 


Zustrom entsprechend einer Staupunktstromung (w=ax, v=—a y) erfolgt. Die 
Reibungsspannung t, am Umfang ist iiberall von gleichem Betrag und folgt 
nach der Theorie zu : 

T =1,15 Vav 


(ey SE = 
1 "(e/2) W? as 


wo W die Umfangsgeschwindigkeit des Zylinders ist. Bei unserem Versuch 
(Figur) floss die Luft aus einer Parallel-Ringdiise, die von der Zylinderoberflache 
den Abstand h hatte. Fiir diesen etwas veranderten Zustrom ergibt sich 
V 
= i ae 
a=O75 h? 
wo V die (radiale) Geschwindigkeit in der Miindung der Ringdiise bedeutet. 
Fiihrt man noch das Verhaltnis B =V/W ein, so ergibt sich 


iirese 
C,=1,015 ae VB, 


_ welcher Wert in guter Ubereinstimmung mit den Messungen ist, falls V geniigend - 


gross ist. 4 
Bei kleinem f tritt Ubergang zur Turbulenz ein. Man kann vermuten, dass 


durch Zentrifugalkrafte zuerst unstabile, aber geordnete Wirbel mit Achsen 
ungefahr in peripherer Richtung entstehen, wie sie von TayLor und GORTLER 
naher betrachtet wurden. — Das Kriterium 

We oO 


v tor ial 


wo # die Impulsmangeldicke und R der Zylinderradius ist, wird mit # ~ 2 erfiillt, 
was plausibel erscheint. Der Zerfall dieser G6rtler-Wirbel wiirde dann zur eigent- 


lichen Turbulenz fiihren. 


Beispiele zur Theorie der Fliigelgitter. Von L. MEYER, Ziirich*). 

Zur Priifung der am Institut fiir Aerodynamik der ETH entwickelten, aut 
reibungsfreier Strémung basierenden Berechnungsmethode fiir Schaufelprofile 
von axialen Turbomaschinen wurden Messungen an Schaufelgittern im elektro- 
lytischen Trog durchgefihrt. 

Nachdem eine geniigende Messgenauigkeit erreicht war®), wurden mehrere 
verschiedenartige Profile untersucht. Die hier wiedergegebene Figur zeigt zum 
Beispiel die gerechneten und gemessenen Geschwindigkeiten an einem 5% dicken 
Profil fiir Axialverdichter bei stossfreiem Eintritt (B,= 26,2°) sowie bei geander- 
ter Anstrémrichtung. Es konnte eine befriedigende Ubereinstimmung zwischen 


Rechnung und Messung festgestellt werden’). 


4) Institut fiir Aerodynamik der ETH. 4 
5) L. Meyer, Zur Theorie der Fligelgitter, Z. angew. Math. Phys. 8, Fasc. 4 (1957). 
6) Ein ausfithrlicher Bericht wird in den «Mitteilungen aus dem Institut fiir Aerodynamik an 


der ETH», herausgegeben von Prof. Dr. J. ACKERET, erscheinen. 
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Profil fiir Axialverdichter. Gerechnete und gemessene Geschwindigkeitsverteilungen bei verschie- 
denen Anstrémwinkeln; die Bezugsgeschwindigkeit V, ist die zur Gitterrichtung normale Kompo- 
nente der Eintrittsgeschwindigkeit. 


Uber die Beliftung von Autostrassen-Tunnel. Von A. HAERTER, Ziirich’). 

Zur geniigenden Verdiinnung der Fahrzeugabgase in der Tunnelluft sind 
Strassentunnel mit Liiftungsanlagen auszustatten. Beim System der Querliiftung 
wird aus einem parallel zum Verkehrsraum laufenden Zuluftkanal Frischluft in 
den Verkehrsraum eingefiihrt und die Tunnelluft wieder in einen Abluftkanal 
abgesogen. 


*) Institut fiir Aerodynamik der ETH. 
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Die Druck- und Geschwindigkeitsverlaufe in den Luftkandlen berechnen sich 
nach untenstehendem Gleichungssystem: 


an, if dw? 
dE = FA 16) O5,— io de ? 


dé ease F/L Vi m,t+Bo2F my, 
wobei 
ong ot Mas. bib heen 
gE ia? Og, = tea ? Nag = (e/2) uz, ? 


4 L/D Rohrreibungsfaktor, f Schlitzéffnung pro d&, F Kanalquerschnitt, ¢ Um- 
lenk- und Reibungsverlustfaktor fiir Schlitz, K, B, von Grenzschicht abhangig, 
K, von v,/u, abhangig (v Einstr6mgeschwindigkeit), , Abluftkanal, , Verkehrs- 
raum, , Zuluftkanal. 

Die Randbedingungen lauten €=0, w=1; €=1, w=0. 

Fir den Fall konstanter Schlitz6ffnungen und konstanten Kanalquerschnitts 
ist die Mengenverteilung langs des Kanals einerseits von 4 L/D und anderseits von 


HVI+E 
Bir s2 


das heisst vom Druckniveau des Kanals gegen aussen abhangig. Wird eine einiger- 
massen gleichmassige Verteilung verlangt, so muss der mittlere Druckunterschied 
gegen aussen im allgemeinen gross sein. 

Durch Einstellen der Schlitze kann natiirlich eine gleichmassige Verteilung 
auch fiir kleine Druckunterschiede gegen aussen erzielt werden, wie es bei der 
Tunnelventilation zwecks Leistungsersparnis gemacht wird. Beobachtungen am 
Schelde-Tunnel in Antwerpen zeigten aber, dass beim Betrieb mit kleinerer Luft- 
menge als jener, fiir welche die Schlitze eingestellt wurden, die Kanale die Luft 
nicht mehr gleichmassig verteilten. 

Dafiir ergeben sich nun zwei Erklarungen: 

Eine Ursache liegt in der Anderung der Reynolds-Zahl entsprechend der 
Mengenanderung, sie bringt aber bei Verhaltnissen, wie sie bei der Tunnelventi- 
lation vorliegen, nur kleine Abweichungen. 

Die Hauptursache ist die, dass der Druck im Verkehrsraum nicht immer 
konstant ist, wie fiir die Schlitzeinstellung angenommen wurde. Atmospharische 
und bei Tunnelsteigungen auch thermostatische Einfliisse bewirken Druckgra- 
dienten im Verkehrsraum,; und der Luftwiderstand der durchfahrenden Fahr- 
zeuge, der sich im Tunnel in einem Druckunterschied vor und hinter dem Wagen 
dussert, ergibt merkliche Druckveranderungen. Am Modellversuch wurde durch 
Kraftwaigung und Druckabfallmessungen der Einfluss der Tunnelwande und des 
Fahrzeugabstandes auf den Luftwiderstand auch gemessen. Diese Druckande- 
rungen k6onnen sich nun speziell bei kleinen Luftmengen auswirken, die nur noch 
unter geringem Uber- bzw. Unterdruck gegen den Verkehrsraum stehen. 

Um den Einfluss der Druckgradienten allgemein tiberblicken zu k6nnen, 
wurde im Bereich der méglichen Kanalabmessungen (durch A L/D charakterisiert) 
und der méglichen Druckunterschiede ,, zwischen Kanalanfang und -ende fiir 
verschiedene Unter- bzw. Uberdriicke 2,,, am Kanalende der maximale Unter- 
schied im Vergleich zur mittleren Menge 4Ao,, berechnet. Die umfangreichen 
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Rechnungen konnten auf der Contraves-Analogie-Rechenmaschine durchgefuhrt 
werden. Ein Beispiel der aufgestellten Kurvenblatter gibt unsere Figur. Mit 
diesen Berechnungen wird es somit moglich sein, Beliiftungskanale optimal aus- 
zulegen. 
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Druckgradienten im Verkehrsraum 7,  (Druckunterschied tiber gesamte Kanallange, bezogen auf 
den Eintrittsstaudruck im Zuluftkanal) ergeben Abweichungen Ao, von der gleichmassigen Vertei- 
lung, je nach Kanaliiberdruck z,, (Uberdruck am Kanalende). 


Uber die Méglichkeit magnetohydrodynamischer Kreisel. Von P. MaE- 
DER, Ziirich®), 


Eine ausfiihrliche Veréffentlichung in der ZAMP erfolgt spater. 


Systematische Herleitung logischer Schaltungen. Von H. GLArtrt, 
Ziirich®). 
In vereinfachter Darstellung besteht das vorgeschlagene Verfahren aus fol- 
genden Schritten: 
1. Systematische Komposition von Elementar-Schaltkreisen. 
2. Klassierung der erhaltenen Schaltkreise zwecks Erstellung eines Katalogs der 
mit Kompositionen aus Rk < R, Grundschaltungen darstellbaren Funktionen. 


8) Zurzeit Institut fiir Aerodynamik der ETH. 
®) IBM-Forschungslaboratorium. 
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3. Abtastung der in Frage kommenden Teile des Katalogs entsprechend der zu 
realisierenden Schaltfunktion; Versuch, die im Katalog enthaltenen Funktio- 
nen mit der zu realisierenden in Deckung zu bringen, notigenfalls durch 
Gleichsetzung oder Fixierung von iiberzahligen Variabeln. 


Ausgehend von der Tatsache, dass alle Vorgange in Schaltkreisen, die kapa- 
zitive und induktive Elemente enthalten, einzeln durch Gleichstromschaltkreise 
dargestellt werden kénnen, werden 6 Elementarschaltkreise, bestehend aus je 
einem nichtlinearen Element (Diode, pnp-Transistor, npn-Transistor oder Hoch- 
vakuumtriode) in Verbindung mit einem Ohmschen Widerstand, eingefiihrt. 

Alle Anschliisse werden als variable Eingadnge aufgefasst. Spannungen werden 
nur qualitativ angegeben. 

Die logischen Funktionen von Kompositionen lassen sich in Boolscher Algebra 
durch einfache Substitution herleiten. Ein geeigneter Bezeichnungscode liefert 
Angaben iiber Art, Anzahl und Verkniipfung der verwendeten Elementarkreise 
sowie tiber die Impedanzverhidltnisse. 

Zur Klassierung der dusserst allgemeine Funktionen darstellenden Komposi- 
tionen werden Abzahlmethoden herangezogen. Sie erweisen sich deshalb als 
erfolgreich, weil jede Variable ein- und nur einmal sowie nur in positiver oder 
nur in negierter Form auftritt. 

Eine zu realisierende, in Boolscher Algebra vorgegebene Funktion wird nach 
ahnlichen Gesichtspunkten durch «Koordinaten» charakterisiert, wie sie bei der 
Klassierung gelten. 

- Gesetzmassigkeiten des Verfahrens gestatten eine rasche Bezeichnung von 
relativ eng begrenzten Katalogteilen, welche potentielle Lésungen enthalten 
k6nnen. 

Demzufolge gestattet das Verfahren Aussagen der folgenden Art: 

«Zur Realisierung einer bestimmten logischen Funktion bestehen fiir eine vor- 
gegebene Maximalzahl R, von nichtlinearen Elementen folgende und nur folgende 
Méglichkeiten. » 

Dieses vereinfachte Verfahren wird fiir den Fall R,)=2 simuliert. Als erstes 
Beispiel weist X=AB’+ A’B («exclusive or») zwei Loésungen auf, wahrend 
X=A(B+C)+BC eine einzige Lésung hat. Nebenbei sei bemerkt, dass fiir 
X=AB+4A’B’ als einzige Funktion mit zwei Variabeln keine Losung fiir 
eg OXISLICT 1, 

In einer verfeinerten Form des Verfahrens werden Gesetzmassigkeiten des 
Aufbaus logischer Funktionen dazu verwendet, um nur noch begrenzte, durch 
das Problem bestimmte Teile des Katalogs zu entwickeln. Dies erméglicht die 
Behandlung komplizierterer Probleme. 

Als Besonderheiten der Methode sind zu erwahnen: Vollstandigkeit (Garantie 
fiir Lésung mit geringstem Aufwand) sowie die Méglichkeit der maschinellen 
Durchfiihrung (hinreichend exakte Formulierung). 


Ein einfacher Thermostat mit Transistoren. Von E. BaLpINcER und 
A. MatrEr, Basel?®). 

Die Temperaturabhangigkeit der Transistoreigenschaften stort in manchen 
Anwendungen, speziell wenn hohe Anforderungen an Genauigkeit und Konstanz 
gestellt werden. Wohl ist durch geeignete Schaltungen eine weitgehende Kom- 
pensation von Temperatureffekten méglich. Haufig wird aber die Lésung eines 


10) Physikalische Anstalt der Universitat, Abteilung fiir angewandte Physik. 
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Schaltungsproblems wesentlich erleichtert oder gar erst erméglicht, wenn ein 
Teil der Schaltung in einem Thermostaten untergebracht werden kann. 

Wir haben einen einfachen Thermostaten entwickelt, der in Betriebsspannung 
und Betriebssicherheit einer Transistorschaltung angepasst ist (Figur 1). Als 


0C76 0C76 OCI6 


Figur 1 
* Schaltschema des Transistor-Thermostaten. Mit R, = 150 Q ergab sich eine Gleichgewichtstempe- 
ratur von 46°C. R, = 10Q. 


Temperaturindikator dient der Transistor T, mit festgehaltener Basisspannung"’). 
Bei fallender Temperatur sinkt sein Kollektorstrom. Dadurch wird der Heizstrom, 
der durch 7; gesteuert wird, vergréssert. T, reduziert, als Emitterfolgestufe, die 
Belastung von 7T,. Mit Hilfe der Zener-Diode gelingt es, die Gleichgewichtstem- 
peratur des Thermostaten von der Speisespannung unabhingig zu machen: R, 
wird gegentiber dem Innenwiderstand der Zener-Diode derart justiert, dass sich 
die Spannungsschwankungen an den Punkten A und B¥2) in ihrer Auswirkung 
auf den Heizstrom so weit kompensieren, dass die gesamte Heizleistung (bei 
gleichbleibender Temperatur) sich nicht andert. 

Die Versuchsausfiihrung des Thermostaten besteht aus einem 3 mm starken, 
beidseitig geschlossenen Messingzylinder von 5 cm Durchmesser und 5 cm Hohe. 
Im Innern ist die Regelschaltung (Figur 1) ausser dem Heizwiderstand von 10 Q 
untergebracht, der auf dem Messingzylinder aufgewickelt ist. Der nutzbare Innen- 
raum betragt 60 cm®, was jedoch keine obere Grenze darstellt. Die Warmeabgabe 
nach aussen ist durch eine allseitig 2 cm dicke Schicht von Polystyrolschaum auf 
75 mW/°C reduziert. 

Der Thermostat wurde mittels R, auf eine Gleichgewichtstemperatur von 46°C 
eingestellt. Der Temperatur-Stabilisierungsfaktor betrug 50 und der Spannungs- 
bereich 8,4 V +12 V. Die Kurzzeitkonstanz iiber 4 Stunden war besser als die 
Messgenauigkeit von + 0,015°C. Uber die Langzeitkonstanz orientiert Figur 2. 
Bei der Beurteilung des Driftes von 0,2°C im Dauerversuch von einem halben 
Jahr ist zu beachten, dass ausser der Alterung der Transistoren und der Zener- 
Diode auch eine Alterung der verwendeten Kohlewiderstande (Resista, 1/2 Watt) 
eine Rolle spielen kann. Die Streuung der Messpunkte diirfte vom Einfluss der 
Umgebungstemperatur auf den einfachen Referenzthermostaten herriihren, Dieser 
wies einen schlechten Temperatur-Stabilisierungsfaktor (rund 20) auf, seine 
Langzeitkonstanz war jedoch durch die Verwendung eines Hg-Kontaktthermo- 

3) 


. E. BALDINGER und W. Czayja, Z. angew. Math. Phys. 9, 22ff. (1958). 
) 


Unter Umstanden miissen die Anschliisse bei_A und B vertauscht werden, 
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meters gesichert. Obwohl die Messresultate, gemiiss der festgestellten Korrelation 
mit der Zimmertemperatur, korrigiert wurden, bleiben infolge der Tragheit der 
Anordnung Restfehler bei schnellen Anderungen der Zimmertemperatur. 
Eine Verbesserung der Temperaturstabilisierung des Transistor-Thermostaten 
_ kann durch Erhéhung der Verstaérkung der Regelschaltung erreicht werden. Ein- 
fache Méglichkeiten sind: 


a) Ersatz von TJ, durch eine Verstarkerstufe mit einem npn-Transistor; 
b) Einfiigung eines kleinen gemeinsamen Emitterwiderstandes fiir T, und TJ; zur 
Erzielung positiver Riickkopplung. 


Nov. Dez. | Jan. Feb. Marz April 


Figur 2 
Gang ATg der Gleichgewichtstemperatur des Transistor -Thermostaten. 


Versuche dieser Art sollen in Verbindung mit einem besseren Referenzthermo- 
staten unternommen werden. 

Fiir die finanzielle Unterstiitzung dieser Untersuchungen sind wir der Stiftung 
Hasler-Werke zu Dank verpflichtet. 


Ein stabilisiertes Netzgerat hoher Giite mit Transistoren. Von E. Bat- 
DINGER und W. Czaja, Basel*). 

Bereits friiher haben wir iiber ein stabilisiertes Netzgerat mit Transistoren 
berichtet, das einen Stabilisierungsfaktor von rund 1000 und einen Innenwider- 
stand von etwa 1 mQ aufwies!4). _ 

Eine Verbesserung dieser Daten ist dann mdglich, wenn eine Verbesserung der 
Referenzspannung gelingt}*). 

Im folgenden wird iiber ein stabilisiertes Netzgerat mit Transistoren berichtet, 
dessen Stabilisierungsfaktor bei einer Ausgangsspannung von 6 V etwa zwischen 
0,5 bis 1,0-104 (ohne Kompensation) liegt und das einen Innenwiderstand von 
etwa 0,2 mQ aufweist. Die Ausgangsspannung ist variabel zwischen 1 und 17 V, 
und die Ausgangsleistung betragt maximal 15 W, wobei ein Laststrom von 2,5 A 
nicht iiberschritten werden darf. Ferner ist das Gerat kurzschluBsicher [die Daten 
sind gleich wie in 15)]. Eine Ausgangsspannung von beispielsweise 6 V ist wahrend 
1 Stunde besser als --0,1 mV konstant (Figur 1). Die Rippelspannung am Ausgang 
betragt weniger als 0,1 mV. Das gesamte stabilisierte Netzgerdt weist etwa die 
Dimensionen eines 1-kW-«Soerensen »-Stabilisators auf. 


18) Physikalische Anstalt der Universitat, Abteilung fiir angewandte Physik. 


14) —, BALDINGER und W. Czayja, Z. angew. Math. Phys. 8, 324 (1957). 
15) —. BALDINGER und W. Czaja, Z. angew. Math. Phys. 9, 1-25 (1958). 
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Figur 2 zeigt das Blockdiagramm des Gerates. Als Thermostat wurde ein 
Transistor-Thermostat*) mit einer Olfiillung verwendet. Zur Vermeidung un- 
kontrollierbarer Temperaturdifferenzen im Apparat selber ist ein Geblase zur 
Kihlung eingebaut. 


© Ausgang 


Netz oO 
Ausgangspannung = 6 Y= 
0,2mV ‘ 
Vorstabi- Referenz — i 
a | | lisierung spannung : 
ql Y; Ys Y; We 
— Stunden —= See 
Figur 1 Figur 2 
Vorlaufige Messung der Ausgangsspannung. Blockschema des stabilisierten Netzgerates. 
Vergleich mit 6-V-Bleiakkumulator, © 9 ++ deutet die im Thermostat untergebrach- 


ten Schaltungen an. 


ae IY AUSGANG 
25004F 
ch : 
ohh a aoe | 
P| 
ree | 
14V= O) O) © 0C76 Philips 
Be bie ee ne -Q@- OCI6) Fae 
| I & a a A @ 2N277 Delco 
® Zs Jntermetall 
® zé ‘: 
ce ®© Fo7 * 
33 
® Figur 3 


Schema des stabilisierten Netzgerates. 


Die Transistoren und Widerstande der Schaltung in Figur 3 sind nicht aus- — 
gesucht. Es wurde lediglich der symmetrische Verstarker durch passende Wahl 
der Spannungsteiler bei der Arbeitstemperatur des Thermostaten (etwa 40°C) 
symmetriert. In der Schaltung wurden ausschliesslich Drahtwiderstande verwen- 
det und Elektrolytkondensatoren nur dort, wo eine Verdnderung tolerierbar ist. 


Diese Untersuchungen wurden aus Mitteln der Stiftung Hasler-Werke finan- 
ziert, wofiir wir unseren Dank aussprechen. 


16) E. BALDINGER und A. Mater, siehe voranstehenden Bericht. 
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‘Ein supraleitender Gleichrichter und Verstarker. Von J. L. OLSEN, 
= Ziirich?’). 

Viele Experimente bei der Temperatur des fliissigen Heliums erfordern relativ 
_grosse Gleichstréme bei sehr niedrigen Spannungen. Um solche Stréme zu liefern, 
ohne dass zu grosse Warmemengen in den Kryostaten eingefiihrt werden, ware 
_ ein in Helium arbeitender Transformator mit einem passenden Gleichrichter 
zweckmassig. Jedoch fehlte bisher ein einfacher, bei tiefen Temperaturen arbei- 
tender Gleichrichter. 

Wir haben jetzt die in Figur 1 gezeigte Schaltung zur Lésung dieses Problems 
untersucht’). Der Sekundarstrom J, des Transformators T wird durch die beiden 


l, 
H, 
H, 
Figur 1 
Schaltbild. 


- aus supraleitendem Material gewickelten Spulen S, und S, gefiihrt. Ubersteigt J, 
einen gewissen Grenzwert J,, so wird in S, und S, wahrend Teilen der Periode die 
Supraleitung zerstért und ein Anstieg von J, iiber J, erschwert. Durch Anlegen 
eines axialen Magnetfeldes auf S, und S, wird der Absolutwert von J, von der 
Stromrichtung abhangig, und mit passender Wahl der Richtung des von den 
Feldspulen erzeugten Feldes wird eine Gleichspannung zwischen den Punkten A 
und C entstehen. 

In einer kiirzlich untersuchten experimentellen Ausfiihrung dieser Schaltung 
hatte der Transformator T 1000 Primarwicklungen und 4 Sekundarwicklungen. 
Der ganze Sekundarkreis wurde aus Létzinndraht von 0,6 mm Durchmesser 
gemacht, da dieses Material etwa die supraleitenden Eigenschaften von Blei, aber 
viel giinstigere mechanische Eigenschaften besitzt. S; und S, bestanden je aus 


Vp 
$f \ 
mV 

i _ 

03 F 
Dac ile, 

0 1 1 ft i = 

0 05 A 10 01 10 A 10 


Figur 3 


Figur 2 
V, in Funktion von Jy fiir 4c = 110 mA. 


V in Funktion von [4c fir Iy = 0,5 A. 


17) Institut fiir kalorische Apparate und Kaltetechnik der ETH. 
18) Siehe auch J. L. Orsen, Rev. Sci. Inst. 1958 (im Druck). 


294 Varia — Miscellaneous — Divers ZAMP 


10 dicht aneinander plazierten Windungen von 2 mm Innendurchmesser. H; und 
H, hatten 40 Windungen pro Zentimeter. 

In den Figuren 2 und 3 zeigen wir die Ausgangsspannung V, zwischen A und C; 
wobei wir in Figur 2 den Feldstrom I, konstant auf 0,5 A und in Figur 3 den 
Primdrstrom I4¢ des Transformators konstant auf 110 mA halten. 

Es ist klar, dass dieses Gerat ausser als Gleichstromgenerator auch als Ver- 
stiirker verwendet werden kann. Mit supraleitenden Feldwicklungen, einer Ein- 
gangsimpedanz von 2-10-*Q und einem Belastungswiderstand von 2-10~¢ Q 
erzielten wir einen Stromverstarkungsfaktor von 10 und eine Spannungsverstar- 
kung von 200. : 

Diese Arbeit wurde aus einem Arbeitsbeschaffungskredit des Bundes finanziell 
unterstiitzt. 


Ein nichtlinearer Verstarker zur Verkiirzung von Impulsflanken. Von 
W. Beuscu, D.MaEpER?®) und R. MULLER, Ziirich”®). 

Bei der Spektroskopie von y-Quanten mit anorganischen Szintillatoren liefert 
die Anode des Photo-Elektronenvervielfachers einen mit der Zeitkonstante 
t = 0,3 us (fiir NaJ) exponentiell abklingenden Strom. Fiir Koinzidenzmessungen 
hoher Zeitauflésung muss ein kleiner Bruchteil dieses Stromimpulses ein Signal 
auslésen, das heisst im wesentlichen eine Gitterkapazitat aufladen. Der hier be- 
schriebene Verstarker steuert eine Endstufe innert 4 ns aus, wenn die Photoréhre 
impulsmassig mindestens 5 yA liefert. Er gestattet hohe Koinzidenzauflésung 
auch ohne Elektronenvervielfacher mit extremer Verstarkung. 

Die Anstiegszeit des ganzen Verstarkers ist wesentlich héher als die genannten 
4 ns?1). Weil der Verstarker in jedem Fall iibersteuert wird, bestimmt die Zeit, 
welche er fiir die ersten 10% des Anstiegs bendtigt, die minimale Dauer des 
Koinzidenzsignals. Der Verstarker besteht aus 8 Rodhren, welche je 1,3- bzw. 
8fach verstarken. Die Anstiegszeit fiir eine einzelne Anode ist ungefahr RC =5 ns. 
Waren die Arbeitswiderstande rein ohmisch, so wiirde der 8stufige Verstarker fiir 
die ersten 10% des Anstiegs 25 ns, fiir die nachsten 80% weitere 30 ns benGtigen, 
ware also fiir unsere Zwecke unbrauchbar; vor allem, weil die Zeit fiir die ersten 
10% noch stark von der Impulsamplitude abhangt. 

Die Induktivitaten L, und die Trimmerkondensatoren bilden mit den Anoden- 
und Gitterkapazitaten ein doppeltes a-Glied einer Laufzeitkette. Die Kapazitaten 
werden teilweise kompensiert, und es resultiert eine von der Impulsamplitude 
unabhangige Laufzeit im Verstarker. Die Wirkungsweise des Kopplungsgliedes 
kann so veranschaulicht werden: Die beiden Spulen von L, bilden einen Trans- 
formator, der einen zum Beispiel negativen Impuls unverzégert mit positivem 
Vorzeichen iibertragt. Dazu wird der ganze [mpuls mit unverandertem Vorzeichen, 
aber mit der Verzégerung des doppelten a-Gliedes iiberlagert. Das resultierende 
Uberschwingen ist in Figur 2, oben, sichtbar, desgleichen die Tatsache, dass die 
ersten 10% des Anstiegs jetzt am raschesten durchlaufen werden. 

Berechnungen haben ergeben, dass der Verlauf der ersten 10% des Anstiegs 
ziemlich kritisch von der Kopplung abhangt und dass k=0,7 der giinstigste 
Wert ist, was durch das Experiment bestatigt wurde. 


19) Zur Zeit am ORNL., Oak Ridge, Tennessee. 
20) Physikalisches Institut der ETH. 


’ *1) Ein Kettenverstarker mit 3-10 Réhren kénnte zum Beispiel bei 1000facher Verstarkung 
diese Anstiegszeit besitzen. 
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Bei der Figur 2 wurden auf den Eingang des Verstarkers exponentiell mit 
t = 0,3 ps ansteigende Impulse der Amplitude 0,2 V gegeben. Eine Vergrésserung 
der Amplitude um einen Faktor 10 veranderte die Form des Impulses nur un- 
wesentlich, er erschien um 5 ns friiher. 


30V 


Figur 2 
Oszillogramme der Impulse am Ausgang. Oben: Impulsform mit «Uberschwingen» nach der 
4. Doppelstufe. Das Uberschwingen wird durch die nachfolgende Seriediode nicht ubertragen. 
Unten: Impulsform nach der Endstufe. Zur Kontrolle der Zeitbasis wurde die Lange des Recht- 
eckimpulses mit einem kurzgeschlossenen Kabel auf 20 ns begrenzt. 


Un amplificateur pour spectrophotométre. Par R. DEessoutavy, J. Du- 
RAND et E. LUscuer, Lausanne et Herisau??). 

Comme décrit?’), le spectrophotométre Metrohm utilise un faisceau lumineux 
interrompu périodiquement et mesuré par une photocellule a vide. 

Il s’agit d’avoir un amplificateur 4 courant alternatif dont le gain peut étre 
ajusté en 8 pas de 10 db. La précision requise est de +0,5% dans toute la gamme 
de lecture. La sensibilité doit étre limitée par le bruit de la résistance d’entrée R. 


Rapport signal | bruit 


Le phototube peut étre assimilé a une source de courant J traversant la résis- 
tance de charge FR. Le rapport signal/bruit est: 


U, RI VR F 
Up V4kTAfR Vat — 


Il y a donc intérét & avoir R aussi élevé que possible. Une limite de R est 
toutefois fixée par l’effet des capacités parasites, effet d’autant plus faible que 
la fréquence est plus basse. Nous avons choisi R = 400 MQ et #/ =10 Hz. Avec ces 
valeurs, le bruit de la lampe d’entrée est négligeable (R équivalant ~ 1 MQ)?4). 


Largeur de bande 


A une fréquence aussi basse, de simples filtres RC, passe-bas et passe-haut, 
permettent d’obtenir une largeur de bande de l’ordre de 12 Hz. 

L’amplitude du signal étant généralement beaucoup plus grande que celle du 
bruit, il n’y a pratiquement aucun avantage a utiliser un redresseur en phase?*). 
Laboratoire Metrohm. 

angew. Math. Phys. 8, 223 


) 

WE, (1957). 
24) J. DarBert, Le Vide 52, 200 (1954). 

) ng 


5 
R. A. Smirn, Proc. Electr. Eng. 98, 43 (1951). 
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Sensibilité 


Avec l’amplification maximum, la déviation totale de l’instrument de mesure 
est donnée par un signal de 1-10-12 A. Ce courant correspond a une intensité de 
lumiére de 0,045 ulm (avec un phototube 90 AV Philips et une lampe de tempé- 

-rature de couleur de 2700°K). 

En l’absence de signal et pour l’amplification maximum, le bruit thermo- 
dynamique de la résistance d’entrée provoque une déviation de 3 a 4% de la 
déviation totale. En présence d’un signal suffisant, le bruit ne provoque plus de 
déviation, mais seulement des fluctuations de la valeur lue. 


Montage 


Le redresseur utilise 2 diodes au silicium présentant une grande résistance 

inverse, méme a une température élevée. I] posséde une caractéristique de redres- 

sement extrémement linéaire (écart de linéarité de +0,3% de la valeur maximum) 
car il est alimenté par une source a haute impédance (pentode sous-alimentée). 


Préamplificateur ; 
aniveau reglable  Amplificateur Galvanometre 
de 10en 10db ajustable Redresseur a index lumineux 


Entree des 
phototubes _ 


Alimentation 
stabilisee 


Conclusions 


Cet amplificateur permet de mesurer avec précision des flux lumineux dans 
une gamme d’intensité variant de 1 4 100000. Sa sensibilité est a la limite de ce 


que peut donner un appareil industriel avec cellule a vide. 


Klystron stabilisé en phase. Par P. Kartascuorr, Neuchatel”). : 
Dans une communication antérieure nous avons décrit un spectrographe a 
haute résolution utilisant un klystron stabilisé en phase sur une harmonique 


élevée de l’oscillation d’un quartz?’). 

26) Laboratoire suisse de Recherches horlogeres. 

2 J. Bonanomt, J. DE Prins, J. HERRMANN et P. Kartascuorr, Helv. Phys. Acta 30, 290 
(1957). 


298 Varia — Miscellaneous — Divers ZAMP 


La fréquence du quartz est balayée par une réactance variable (figure 1). La 
self L, varie grace a la saturation du noyau en ferrite. Une polarisation appropriée 
permet d’obtenir une bonne linéarité du balayage. 


cam -ooles» ferpoxcube Lb 3 
Polarisation = Type D 14-8 


s 4 Excitation 
WW 


Oscillateur 4 quartz balayé en fréquence par une self variable. 


Alimentation 
Klystron 


Batterie 
B00V 


Figure 1 


Klystron X13 


Doubleur 
a diode 
Q= 


Sortie 23870 MHZ 


Tension 
de correction 


Oscillateur 
auxiliaire 


Figure 2 
Circuit de stabilisation en phase d’un oscillateur klystron. 


Pour que le domaine de maintien de la synchronisation du klystron soit plus 
grand que celui du balayage maximum il faut que le gain de la boucle d’asservis- 
sement soit suffisant. Les propriétés du filtre de la tension de correction (figure 2) 
permettent un gain plus élevé que dans le circuit original de PETER et STRAND- 
BERG?8), 

Nous avons obtenu sans difficulté un domaine de synchronisation de 7 MHz 
autour de la fréquence de 24000 MHz. Le balayage était de quelques centaines 
de kHz. 


28) M. Perer and M. W. P. SrraNnpBERG, Proc. Inst. Radio Engrs. 43, 869 (1955). 
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La fréquence auxiliaire était de 150 kHz mais elle peut étre choisie entre 
100 kHz et plusieurs MHz, la limite supérieure étant donnée par Vamplificateur 
et le discriminateur de phase. Pour des valeurs inférieures A 100 kHz le domaine 
de synchronisation diminue. 

La tension de correction est appliquée au réflecteur du klystron a travers une 

‘batterie qui fournit la tension fixe du réflecteur. C’est dans ce sens que ce 
systeme différe de celui réalisé par Davis®*) ot le discriminateur de phase et le 
filtre se trouvent au potentiel élevé du réflecteur. 


Summary 


In an earlier communication on High Resolution Microwave Spectroscopy 2’) 
we described a spectrograph using a phase-stabilized klystron oscillator as a local 
oscillator in the detection system. The stabilizing system has a ‘holding range’ 
of several megacycles which permits the sweeping of the reference crystal oscil- 
lator. 


Multiplication de fréquence et stabilité en phase. Par P. KARTASCHOFF, 
Neuchatel?$). 

Les oscillateurs a quartz 100 kHz fournissent une fréquence qui est trés stable 
a long terme. Quand on multiplie la fréquence pour arriver a des fréquences de 
10000 MHz ou plus on observe que l’instabilité de la phase élargit la raie spectrale 
du signal. 

Par contre il a été démontré que des oscillateurs a quartz relativement simples 
qui travaillent sur des fréquences de 5 MHz ou plus donnent un signal pur et 
stable 4 mieux que 10~1° a court terme*®). Ces performances sont réalisées avec 
des quartz de télécommunication ordinaires (coupe AT). Si le quartz est thermo- 
statisé la dérive lente de la fréquence est inférieure 4 10~8 par jour, mais toujours 
plus grande que celle d’un oscillateur étalon 100 kHz. oe 

S’il est possible de maintenir la fréquence d’un oscillateur HF sur un multiple 
exact de 100 kHz on peut obtenir un signal de fréquence plus stable et dont la 
raie spectrale est pure. 

Ceci est réalisé par un systéme d’asservissement de la phase (figure 1) qui est 
basé sur le principe IGO#*1)?), 


Oscillateur Oscillateur a quartz 
étalon 100 KHz haute fréquence 


| Formation Diser. 
impulsions de phase 


Figure 1 
Circuit de stabilisation par impulsions d’un oscillateur (TGO). 


Sortie 


Correction 
fréquence 


29) B. F. Davis, External Publication No. 380 (Jet Propulsion Laboratory, California Inst. of 
Technology, 1957). 

30) M. W. P. SrranpBERG, Letter Proc. IRE 44, 696 (1956). 

31) —. G. HucENuortz, Philips techn. Rundschau 13, 11 (1952). 

32) G. Sarmet, Proc. IRE 44, 1582 (1956). 
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La phase de l’oscillation HF est comparée a des impulsions 100 kHz etal 
signal de correction passe a travers un filtre intégrateur (figure 2) qui a une 
constante de temps élevée (10-20s) afin d’éliminer toutes fluctuations bréves. 
Pour la méme raison le domaine de maintien de l’asservissement doit étre aussi 
petit que possible. Dans nos essais il était de quelques 10~* ce qui suffisait pour 
produire un signal propre et pour maintenir la fréquence asservie pendant des 
heures. 


100mH 10MQ 


é IMQ ; 
Entree Sortie 


Figure 2 


Filtre intégrateur qui élimine l’effet des fluctuations breves de la phase. 


Summary 


A frequency multiplying system which delivers an exact multiple of the 
frequency of a 100 kc/s standard in the microwave region has been developed. 
High spectral purity of the signal has been obtained by servo-controlling the 
phase of a high frequency quartz crystal oscillator through an integrating filter 
with long time-constant. 


Verwendung einer Differentialphotoréhre in einem hochstabilisieren- 
den Magnetstromregulator. Von D. BRINKMANN, E. BRuN und H. H. Staves, 
Ziirich #*), 

Fir die Beobachtung von Kerninduktionssignalen der Breite von etwa 10 mGs 
bei einem Magnetfeld von rund 10 kGs ist eine Stabilitat (Verhaltnis des Magnet- 
feldes zu seiner Schwankung) von mehr als 10° erforderlich. Der Stabilisierungs- 
faktor (Verhaltnis der relativen Anderungen von Primarspannung und Magnet- 
feld) sollte etwa 10° sein: Eine 1%ige Spannungsschwankung wiirde das Magnet- 
feld um 1 mG dandern. 

Der bisher von uns benutzte Magnetstromstabilisator, bestehend aus Diffe- 
rential- und Gleichstromverstirker (Figur 1), arbeitete mit negativer Riickkopp- 
lung und hatte einen Stabilisierungsfaktor von etwa 3-104 und eine Stabilitat 
von 10°. Durch Verwendung eines Spiegelgalvanometers, das die Magnetstrom- 
schwankungen auf eine Differentialphotoréhre (RCA 5584) tibertragt, wurde der 
Riickkopplungskreis des Stabilisators verbessert und seine Verstarkung herauf- 
gesetzt. Es wurden ein Stabilisierungsfaktor von 105 und eine Stabilitat von 
5-10° erreicht. Ein Magnetfeld von 10 kGs schwankt somit jetzt um etwa 2 mGs. 

Die Aufnahmen der Figur 2 zeigen Dispersionssignale vom Deuterium, an 
denen die Verbesserung des Stabilisators ersichtlich ist. Die oberen Signale 
wurden ohne, die unteren mit Photordhren-Stabilisierung aufgenommen. Bei den 
jeweiligen Aufnahmen rechts ist das konstante Magnetfeld auf die Mitte des 
linearen Teils des Signales eingestellt worden. Das registrierende Schreibgerat 
zeichnet direkt die Schwankungen auf. 


33) Physikalisches Institut der Universitat. 
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Schaltschema des Magnetstromstabilisators. 
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Stabilisation mit Photordhre 


Figur 2 


Dispersionssignale vom Deuterium. 
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Gas Dynamics of Thin Bodies. Von F. I. Franke und E. A. KARPOVICH. 
Ubersetzt aus dem Russischen von M. D. FrreEpMaAN (Interscience Publishers. 
New York und London 1954). 175 S., 33 Fig.; $5.75. 

Der Ubersetzer des russischen Buches iiber die Gasdynamik der diinnen Korper 
erklart im Vorwort, dass die Ubersetzung ins Englische zunachst nur fiir den 
internen Gebrauch im Ames Aeronautical Laboratory der NACA vorgesehen war. 
Die Klarheit und Ausfiihrlichkeit der erklarenden Texte habe ihn jedoch dazu 
bewogen, die vollstandige Ubersetzung zu publizieren. 

Im ersten Kapitel wird betont, dass auch Russland seinen Teil zur Begriindung 
der Theorie der Gasdynamik beigetragen habe, wobei man mit vielen fiir west- 
liche Leser unbekannten Namen vertraut gemacht wird. Anschliessend werden 
die Grundgleichungen hergeleitet und der Begriff des «verzégerten Potentials» 
eingefiihrt, welches sich aus der hier weniger gebrauchlichen Wahl eines mit der 
Stromung mitwandernden Koordinatensystems ergibt. 

Im zweiten Kapitel wird der Geschwindigkeitspotential fiir schlanke Rota- 
tionskérper und ahnliche Korper in Unter- und Uberschallstr6mungen berechnet 
und auch die beschleunigte Bewegung eines solchen Koérpers betrachtet. Hierbei ~ 
ergibt sich das interessante Resultat, dass die Druckkrafte auf einen mit ungefahr 
Schallgeschwindigkeit fliegenden K6érper erst bei Beschleunigungen von tiber 100g 
merklich andern. 

Das dritte Kapitel behandelt im wesentlichen die unendliche und die endliche 
Tragflache in Unter- und Uberschallstro6mungen, wahrend im vierten Kapitel die 
instationare Stro6mung um oszillierende unendliche und endliche Tragflachen, das 
Potential einer pulsierenden Quelle und die Theorie von Hochgeschwindigkeits- 
propellern (mit Unter- und Uberschallgeschwindigkeit der Blattspitzen) bespro- 
chen werden. Hier versagt allerdings die linearisierte Theorie in der Nahe des 
kritischen Propellerradius (Schallgeschwindigkeit), so dass die Resultate wohl als 
theoretisch interessant, praktisch jedoch wenig brauchbar erscheinen. 

Im letzten Kapitel werden die konischen Stré6mungsfelder behandelt. Fiir 
Physiker und theoretisch arbeitende Aerodynamiker diirfte das Buch besonders 
wegen des ungewohnten Einblickes in russische Arbeiten von Interesse sein. Dem 
Konstrukteur mégen die Hinweise fehlen, welche Genauigkeit er von den durch 
die linearisierte Theorie berechneten Resultaten erwarten darf. H.R. VoELLMy 


Praktische Physik. Von F. Kourrauscu, herausgegeben von H. EBERT 
und E. Justi (B. G. Teubner Verlagsgesellschaft, Stuttgart, Bd. 1 1955, Bd. 2 
1956), Bd. 1 646 $., 394 Abb.; DM 36.80, Bd. 2 756 Si, 435 Abb.. 133) Tan 
DM 52.-. 

Fir diese 20. Auflage wurde der altbekannte « Kohlrausch» von den Heraus- 
gebern unter Beiziehung eines erweiterten Mitarbeiterstabes vollstandig neu redi- 
giert. Die Aufteilung des Stoffes auf die beiden Bande ist im wesentlichen gleich . 
geblieben, namlich: 

Band 1: Mechanik, Akustik, Warme, Optik. 
Band 2: Elektrizitat und Magnetismus, Korpuskeln und Quanten; Tabellen zu 
beiden Banden. 

An Neuem gegeniiber friiheren Auflagen*ist besonders hervorzuheben: All- 
gemeine Erweiterung auf allen Gebieten unter Beriicksichtigung der Forschungs- 
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ergebnisse bis Herbst 1954; bessere Ubersichtlichkeit der Tabellen, deren Zahl zu 


-gutem Nutzen bedeutend vermehrt wurde; die Zehnernumerierung der Abschnitte 


tragt viel zum leichteren Auffinden der einzelnen Begriffe bei; das wieder neu 


-aufgenommene Kapitel « Technisches» gibt manchen niitzlichen Hinweis fiir 


praktische Arbeiten im Laboratorium. 


Wie der Zusatz zum Titel (Zum Gebrauch fiir Unterricht, Forschung und 
Technik) bereits andeutet, ist dieses Werk fiir alle in einem physikalischen 


_ Laboratorium Arbeitenden bestimmt. Stellt sich das Problem der Darstellung 


_oder Messung irgendeiner physikalischen Grésse, so kann man sich anhand des 


« Kohlrausch» auf wenigen Seiten iiber die bisher erprobten Methoden rasch 


_orientieren, wobei auch die theoretischen Grundlagen kurz gestreift werden; 
um die Handlichkeit des Werkes nicht zu beeintrachtigen, mussten nattirlich 


allzu detaillierte Angaben vielfach weggelassen werden ebenso wie die Beschrei- 
bung gewisser sehr spezieller Methoden; in diesen Fallen aber wird dem Suchenden 
die Fiille der Literaturzitate (Originalarbeiten und zusammenfassende Darstel- 


- lungen) weiterhelfen. 


Der schwierige Kompromiss, einen Stoff von gewaltigem Ausmass in ein 
wirkliches Handbuch zum Laboratoriumsgebrauch zusammenzudrangen, ist den 


- Herausgebern wiederum gut gelungen. U. Katz 


7 


Technische Schwingungslehre in ihren Grundzitigen. Von E. HUBNER 
(Springer-Verlag, Berlin 1957). 322 S, 208 Abb.; DM 29.40. 

Der Springer-Verlag hat eine Anzahl guter Schwingungsbiicher herausgebracht. 
Man denke etwa an die Technische Schwingungslehre von Hort, an das gleichna- 
mige, leider noch immer fragmentarische Buch von KLOTTER oder an die ver- 
dienstvolle Mesmersche Ubertragung der Mechanical Vibrations von DEN HarRTOoG. 
Wer diese Werke auch nur fliichtig kennt, dem bleibt es unverstandlich, wie der 
Verfasser des vorliegenden Buches (gemass Vorwort) eine Liicke findet zwischen 
den Werken «von hohem wissenschaftlichen Niveau..., die deshalb fiir den 
durchschnittlich mathematisch Geschulten nur schwer zuganglich sind », und den 
iibrigen Lehrbiichern, die «gut verstandlich sind», aber «nicht tief genug in die 
schwingungstechnischen Problemstellungen eindringen». Das vorliegende Buch 
schliesst keine Liicke. 

Dass die Periode einer Schwingung nunmehr « Periodendauer» heisst (S. 1), 
die Abszissenachse kurzweg « Abszisse » (S. 4), dass eine Feder nicht etwa verlan- 
gert, sondern «gelangt» wird (S. 28) und dass eine Beschleunigung «erfolgt » 
(S. 28), sind begriffliche und stilistische Details, die lediglich nebenbei vermerkt 
werden sollen. Schwerer zu iibersehen ist es, wenn etwa (S. 3) ein Winkel im 
Text mit y und in der zugehérigen Abbildung mit » bezeichnet ist, wenn eine 
Summe (S. 10) iiber mehr Summanden erstreckt wird, als in Wirklichkeit vorhan- 
den sind, oder wenn das axiale Tragheitsmoment eines Querschnittes (S. 30) in 
den Rang eines «dquatorialen Tragheitsmomentes » erhoben wird. Dass eine 
Funktion y(t) in der Form y, geschrieben wird (S. 3), ist ja eine leider noch weit- 
verbreitete Gepflogenheit; trotzdem sollte bei einer komplexen Zeitfunktion z(t) 
die reelle Achse (S. 15) nicht mit der Zeitachse verwechselt werden. 

Neben diesen kleineren Mangeln enthalt das Hiibnersche Buch aber auch 
viele grébere Fehler. ; 

So wird zum Beispiel (S. 3) festgestellt: « Winkelfunktionen Se einander 
gleich, wenn ihre Argumente gleich gross sind oder sich um den Winkel 2 x (oder 
ein Vielfaches davon) unterscheiden. » Bei der Behandlung der Schwebung (S220) 
(die er iibrigens auch als « pseudoharmonische Schwingung» bezeichnet) kommt 
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der Autor auf Grund einer Diskussion in der komplexen Ebene zum Ergebnis, 
dass die Frequenz der resultierenden Schwingung stetig verdnderlich sei. Als Vor- 
aussetzung der Entwickelbarkeit einer Funktion in eine Fouriersche Reihe wird 
(S. 21) neben ihrer Endlichkeit ihre Differentiierbarkeit angefiihrt. 

Behandelt man einen gedampften Schwinger mit Hilfe des d’Alembertschen 
Prinzips, so hat man bekanntlich die am Schwinger angreifenden Krafte, namlich 
die Federkraft, die Dampfungskraft und die Tragheitskraft, ins Gleichgewicht zu 
bringen. Manchen Studierenden macht dieser einfache Sachverhalt erfahrungs- 
gemass Miihe. Die Situation wird aber schlechthin hoffnungslos, wenn man, wie 
dies der Verfasser (S. 32, 38) wortlich und (S. 71 usw.) in tibertragenem Sinne 
durch den ganzen Text hindurch tut, die Federkraft als Reaktion der an der 
Feder ausgeiibten Kraft, die Dampfungskraft als Reaktion der « zur Uberwindung 
der Reibungskrafte erforderlichen Kraft» und schliesslich die Tragheitskraft als 
Reaktion einer Massenkraft» definiert, die zur Beschleunigung des Schwingers 
aufzuwenden sei. 

Diese Beispiele diirften ausreichen, um darzulegen, dass das Buch kaum so 
«tief in die Zusammenhange eindringt», wie dies der Autor im Vorwort wortlich 
in Aussicht stellt. Im iibrigen fehlt es auch nicht an Stellen mit einer heiteren 
Note, so zum Beispiel, wenn (S. 133) die Auffassung vertreten wird, dass das 
Verfahren von RAYLEIGH eine Verallgemeinerung des Energietheorems auf 
Systeme mit mehreren Freiheitsgraden darstelle. H.. ZIEGLER 


High Altitude Rocket Research. Von Homer E. NEWELL, jr. (Academic 
Press, Inc., New York 1953). 298 S., 93 Abb.; $7.50. 

Die von den Amerikanern am Ende des Zweiten Weltkrieges in Deutschland 
erbeuteten V 2-Raketen wurden in Amerika teilweise zur Erforschung der héheren 
Atmosphare eingesetzt. Spater konnten vom « Upper Atmosphere- Rocket Research 
Panel» noch die in Amerika entwickelten Héhenraketen Aerobee und Viking 
eingesetzt werden. Die in den Jahren 1946 bis 1952 gesammelten Resultate sind 
im vorliegenden Buch in allgemeinverstaéndlicher Form ver6ffentlicht. 

Die einzelnen Abschnitte des Buches sind nach den folgenden Sachgebieten 
geordnet: Die Rakete als Messk6rper, Mess- und Fernmeldetechnik, Druck, 
Dichte und Temperatur in der hGheren Atmosphare, Winde, Messung der Strah- 
lung der Sonne, Zusammensetzung der oberen Atmosphare, die Ionosphare und 
das erdmagnetische Feld, kosmische Strahlen, Héhenphotographie. 

Die gesammelten Daten erstrecken sich zum Teil bis auf Héhen von 220 km. 
Der Verfasser hat sich verdankenswerterweise die Miihe genommen, nicht nur die 
Messresultate zu verdffentlichen, sondern auch einen kurzen Abriss iiber die 
Messtechnik und die Auswertemethode zu geben. Besonders diese Seite des 
Buches diirfte einen grésseren Leserkreis anziehen, wahrend es den Physikern, 
Meteorologen und Aerodynamikern tiberlassen bleibt, die in Kurven- und Tabel- 
lenform erfassten Resultate weiter zu verarbeiten. H. R. VOELLMY 


Magnetohydrodynamics. Von T. G. CowLinG ([Interscience Tracts on 
Physics and Astronomy, number 4] Interscience Publishers, New York and 
London 1957). 115 S., 17 Abb.; $1.75 / $3.50. 

_ Als G. E. Hae vor einem halben Jahrhundert die Magnetfelder der Sonnen- 
flecken entdeckte, konnte kein Zweifel bestehen, dass viele Erscheinungen der 
Sonnenphysik, zu denen spater die Bewegungen der Protuberanzen und Erup- 
tionen sowie die Struktur der Korona hinzukamen, nur durch koordinierte 
thermodynamische, hydromechanische und magnetische Betrachtungen zu ver- 
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stehen sind. Trotzdem kann erst seit etwa 15 Jahren von einer Magnetohydro- 


_ dynamik gesprochen werden. Wenn auch ihre vornehmlichsten Gefilde Probleme 


der Astrophysik sind, so hat sie in neuerer Zeit doch auch Eingang gefunden in 
der Geo- und Plasmaphysik. Einige wenige der ausschliesslich theoretischen 
Resultate konnten bereits auch durch Modellversuche an bewegten leitenden 
Flissigkeiten in Magnetfeldern verifiziert werden. In diesem Stadium dieser 
jungen Wissenschaft ist das knappe Biichlein von T. G. Cow.rna, der an ihrer 
Entwicklung einen starken Anteil hat, sehr zu begriissen. Der Autor fiihrt im 
ersten Kapitel in die Grundgleichungen der Magnetohydrodynamik ein und 
behandelt auf ihrer Grundlage in den folgenden Kapiteln eine Reihe von speziel- 
len Problemen, wie das magnetostatische Gleichgewicht von Sonnenflecken, die 
Spiralarme der Galaxien, die Sternrotation, die magnetischen Sterne, die magneto- 
hydromechanischen Wellen und ihre Anwendung auf die Theorie der Sonnen- 
flecken, sowie interstellare und terrestrische Magnetfelder. Durch diese Abwand- 
lung des Themas bleibt der Leser dieser anspruchsvollen Monographie bis zum 
Schluss in Spannung. Leider k6nnen stets nur die Wege zu den Lésungen, nicht 
aber diese selbst aufgezeigt werden. Der Grund dafiir liegt einerseits in dem fast 
volligen Fehlen experimenteller Erfahrung in der Magnetohydrodynamik und 
damit im Mangel an Intuition auf diesem Gebiet, besonders aber in den fast 
uniiberwindlichen mathematischen Schwierigkeiten bei der Behandlung nicht- 
linearer Probleme. In zwei Punkten scheint diesem Biichlein eines Pioniers ein 
Erfolg sicher zu sein: in der Klarung der noch vielfach wirren Ansichten auf dem 
Gebiet der Magnetohydrodynamik und in der Weckung von Interesse und Ver- 
standnis ftir dieses dornenvolle Neuland. M. WALDMEIER 


Praktische Mathematik fiir Ingenieure und Physiker. Von R. ZuRMUHL. 
2. Auflage (Springer-Verlag, Berlin 1957). 524 S.,.125 Abb.; DM 28.50. 

Diese bereits in weiten Kreisen bekannte Einfiihrung in die numerischen Ver- 
fahren der praktischen Mathematik befasst sich hauptsaéchlich mit linearen 
Gleichungssystemen, Matrizen-Eigenwertaufgaben, Ausgleichsrechnung, Rand- 
und Eigenwertaufgaben von Differentialgleichungen und anderem mehr. Fiir 
jede der angegebenen Methoden wird die mathematische Fragestellung klar 
herausgearbeitet, wobei zum Verstandnis lediglich die Grundbegriffe der Integral- 
und Differentialrechnung vorausgesetzt werden. Das Hauptgewicht liegt jedoch 
auf der praktischen Anwendung, was mit vielen durchgerechneten Beispielen 
unterstrichen wird. ; 

Fiir die zweite Auflage wurde ausser der Neufassung einiger Kapitel auch ein 


Abschnitt iiber die wichtigsten Priifverfahren der Statistik neu aufgenommen. 
Tu. GINSBURG 


Einfiihrung in die Methoden der Programmierung kaufmannischer 
und wissenschaftlicher Probleme fiir elektronische Rechenanlagen. 
1. Teil: Die Logik der Programmierung. Von BRUNO THURING (Goller Verlag, Baden- 
Baden 1957). 200 S., 100 Abb.; DM 45.-. 

Der vorliegende erste Band behandelt diejenigen elementaren Gesichtspunkte 
der Programmierung digitaler Rechenanlagen, die von den konstruktiven Beson- 
derheiten der jeweiligen Anlage unabhangig sind. Er setzt keine besonderen 


- Vorkenntnisse voraus und ist deshalb als Einfiihrung in die Programmierungs- 


technik fiir National6konomen und Betriebswissenschafter ebenso wie fiir Mathe- 
matiker und Ingenieure geeignet. 
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Im ersten Kapitel werden einige grundlegende Dinge, wie das duale Zahlen- 


system, Buchstaben- und Zahlenverschliisselung und das Rechnen mit festem 


und beweglichem Komma, behandelt. * a 


Im zweiten und dritten Kapitel wird anhand von kaufmannischen (2. Kapitel) 


und wissenschaftlichen Problemen (3. Kapitel) der Aufbau eines Rechenplanes 


auseinandergesetzt, was in erster Linie durch eine ausfihrliche Betrachtung der 
Fluss- oder Ablaufdiagramme geschieht. Als Ubergang zur Befehlsverschliisselung 
wird die algorithmische Schreibweise eines Rechenplanes nach H. RUTISHAUSER 
eingefthrt. : ; 
Das vierte Kapitel befasst sich mit den besonderen Problemen der Ein- und 
Ausgabe, wie sie sich beim Zusammenspiel innerer und ausserer Speicher ergeben. 
Zu bedauern ist, dass friihere Arbeiten iiber die Programmierung (insbesondere 


das grundlegende Werk von J. VON NEUMANN und GOLDSTINE), von denen natiir- 


lich vieles ibernommen wurde, nicht verzeichnet sind. 

Das einfach und klar geschriebene Buch, das sich auch durch eine anspre- 
chende drucktechnische Gestaltung auszeichnet, kann als erste Einfiihrung in 
das Gebiet der Programmierung bestens empfohlen werden. W. NEF 


Prazisionsmessungen von Kapazitaten, Induktivitaten und Zeitkon- 
stanten. Von E. BLEcHScHMIDT. Band 2 (Friedr. Vieweg und Sohn, Braunschweig 
1957). 166 S., 74 Abb.; DM 11.80. 

Anlasslich der Besprechung des ersten Teiles dieser Monographie (Kapazi- 
tatsmessungen) wurde bereits auf das baldige Erscheinen des zweiten erganzenden 
Heftes tiber Induktivitatsmessungen hingewiesen [ZAMP 8, 430 (1957)], das nun 
erschienen ist. Die tibersichtliche Zusammenstellung iiber Induktivitatsnormale 
und tuber die bewahrten Messverfahren zur Untersuchung von Spulen darf Phy- 
sikern und Ingenieuren unbedenklich empfohlen werden. E. BALDINGER 


Momentum Transfer in Fluids. Von W. H. Corcoran, J. B. OpFELL und 
B. H. Sacre (Academic Press Inc., New York 1956). 394 S., 155 Fig.; $9.-. 

Die Verfasser sind von der Absicht ausgegangen, diejenigen Teile der Stré- 
mungslehre zusammenzufassen, die fiir den Chemieingenieur von Interesse sind. 
Daraus erklart sich die etwas eigentiimliche Auswahl und Abgrenzung des 
Stoffes. Deutlich spiegelt sich hier die amerikanische Methode der Ingenieur- 
ausbildung wieder, die von vornherein mehr den Spezialisten anstrebt als die 
europaische. Das vorliegende Buch ist aber nichtsdestoweniger weit davon ent- 
fernt, nur fiir den Chemieingenieur von Nutzen zu sein, sondern es verdient die 
Beachtung aller, die sich mit industriellen Stromungsproblemen abzugeben haben. 

Es werden zunachst die wichtigsten Grundtatsachen behandelt (Bernoullische 
Gleichung usw.). Recht ausfiihrlich werden Erscheinungen der Turbulenz be- 
sprochen. Auch eine Einfiihrung-in die Grenzschichttheorie wird gegeben. In 
einem Anhang werden die Elemente der Tensorrechnung gegeben, und auf dieser 
Grundlage wird in einer mathematisch systematischeren Weise nochmals auf die 
statistische Turbulenztheorie zuriickgegriffen. Das Buch beschrankt sich im 
wesentlichen auf die zusammenfassende Darlegung bekannter Tatsachen, dies 
aber in einer Weise, die dem Studierenden und dem Ingenieur in der Praxis eine 


' wertvolle Hilfe sein kann. Uber gewisse formale Schonheitsfehler, wie etwa dass 


immer noch mit dem spezifischen Gewicht statt mit der Dichte gearbeitet wird, 
kann man wegsehen. W. TRAUPEL 


Vient de paraitre 


— @ 


GARNIR 


Les problemes aux limites de la physique 


mathematique 


Par H. G. Garnir, Chargé de cours & l'Université de Liége (Belgique) 
234 pages avec 5 figures. Prix Fr. 29.— (DM 29.-) 


Mathematische Reihe, Band 23 
«Lehrbiicher und Monographien aus dem Gebiete der exakten Wissenschaften» 


Cet ouvrage concerne l’application des méthodes de l’analyse fonctionnelle 
4 la théorie générale des problémes aux limites de l’équation des ondes et de 
la diffusion. La fonction de GREEN de ces problémes est étudiée minutieuse- 
ment et son réle dans la résolution effective mis en évidence. Ce livre contient 
notamment une justification définitive des méthodes heuristiques basées sur 
les transformations de LapLacr et de FourRIER par lesquelles les physiciens 


recherchent ces fonctions de GREEN. L’auteur fait largement appel a di-' 


verses théories transcendantes comme la théorie actuelle des espaces hilber- 
tiens et la théorie des distributions de L. ScHwart1z. Cependant, l’exposé a été 


rendu self-consistent et présenté de maniére 4 étre compris d’un lecteur pos-- 


sédant seulement une bonne connaissance de l’analyse classique. ° 


' Matiéres: Théorie des espaces fonctionnels hilbertiens — Probléme de DiricH- 
LET-NEUMANN pour l’opérateur métaharmonique — Transformation de 
LAapLace inverse dans la théorie des distributions — Probléme de DrricHLEt- 

NEUMANN pour les opérateurs des ondes et de la diffusion. 
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